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«صححص ‏ مصتهام ه صتطاده وملهعی آهعممصها عم لمتادمه اصعتمن اج امهعماصا هنامرجم تمه نامع مه عمهتتا 
0 )صها0۲مظ1 15 1 ,۰۲3۵۲۵۲۵۲۵ .ممتانه‌تتوه ما ممتاقصطاصصمی ها ۱۱۵۵0 تج وعمتاموع۲ تمه0 نامع 4ص ممهتای رصتقام 
۶ صمتقصوووه ۲2۵10 مطا مااووعنا ,عمصهلهها ماه ممهاه هه همع ۵ فاصمومموصصمم فص عتمتص رامهتناهمه 
۶ ]حمصصووعووه 1۳16812160 2۳0 1۷6مممصم۲م‌م۵ع 10۲ ۵60 2 او ور عمط روع620ع0 07۵ ]ف2 ما کعبن فاعم لهمزعمصناظ 
5 0۱7۷26 ا0ع هه ماو طاهه ۵۶ ممتامعصطمعتعص1 فطع رتقانه‌تانهم ما .عاممصمم‌جصرم عمج صنامع 4ص ممهتاه 
۰ ۲۵0۲ ۱۵۹۵00۲۵۵0 عط) مر ممصهلهه تمه معط رللمزههموه رعماه)صنامع مه ومدیاه طاهها آععصصمع مه ]1۳00۴ 15 
۵۲ 1۳22010۴ ۵۶ مناد فیاممصهآنصته صم ممتاهام0۵6 0680 مه میت معتهطمع۲ ماه ۵۴ 61160 6 تناو وتطا 1 
۵۵۲ م6 2 نطو ۵ ۱۵۵0 ۷/۵۲۵ ۵061 ٩۱۷۷۸۲‏ مط حصم امنهر مه روع)‌یمو تماهصامع فصه ممهتتاه ۳۵ 

۰ ۷۷۲۲ - ۱۷۲۱۱۵۲۲۱۵۱۷۷ 0عصزطامجمی فطه عصتوده عمصفلهه 


۵۹ 2۱0 واهز۱۷]۸)۵ 

م0 1266ناصصله 60 7260مصه 728 8۵061 ۱۷۷/۸۲ عط) حصم فمصتهاهاه فعبله۷ معتقطامع۲ ۵۶ )همع مطا ,تاه فنطا 1 
۱۷۲۵۱۱۳۲۲۸۵۷۷۰ ۵۶ ۳۵01 12)60ع6)ظ1 عمط عصتوتا فاصممممصصمی ممصفلهها ماه 0صنمع مضه رقاطها تمتهه ۵۶ ممتاهتتاهت11 
ممناحعن1 ۵۶ امدمصا فطا مصترفتاصهتن صا دوماو صه)عممصصا فطا ۵۶ فمم) مصتقام. 20طقط ۱۷]2‏ فط صا تع0مصه ۱۷۷ 
عتال 16۵2۲26 12)60اصصه عمط ۴/۵8 معصهاده عع0۱726صامع صرح ممهعداه عط رم 120-986 صا ععصقطه فطل فصه رتصعصصعع‌مصم 
عط و2 ۱۵60 ۷/۵5 مضه ,م۵0 ۸۵۲ 5۱۷۷ معط وطا عم فه ونط 1 .م2۵ صمتامعتا فص المکصنه: ۵۶ ممتادامع ومع معط ما 
0 ۷۷۵۲۵ ۱۱۵061 م۵ صا0ع ۵۶ ممتاهناه۷ مه مم)ه 02 061مصص ۱۷۲۵1۳۲۲۸۵۷۷۱۷۷۲ عطع مه صمتاتجمن رتقک‌صتمه 
من ونهع۷ 10 10۲ م۵0 قفه۳7 ۵۵۱۵0 ممتاهاتاددژه م1 .06مطامصه 1 قاط متاممطماتاج هه 6۲۲۵۲-صه هن اه فظم 
۲ ۵۵۲۱۵0 طا وه 2۵0 ۷۷۵۲۵ وتقع۷ 4 همه 6 طمتطه عم ,2018-2019 60 2009-2010 ۵۶ بقع ۳۶۲01021621 فطا 
۰ ,2018-2019 0 2015-2016 280 2014-2015 10 2009-2010 مصم مه ممتاملتله۷ مه انالوم 
۵ ۱۷۵۵ تاو ما رعع۳۵۵۴0]. .ونقع۷ امن 10۲ امعم دهد میم فاصعمممم‌صمع معصهلهه ماه 0صنام) 
ماه لهتانصا مه میا مه مه ۷۱۵10 لوصو مه الوم متانیهل ۳۱۱ .هه تقصصررمه م24 باه رتتل 10۲ 

1۰ ۱۷۷-۷۷۲۲ (6, ۱۷۲۵۲۲۷۲ ط) ما وتمامصصهتهم 


مکی کز(۱ 20 وااباعع 

مطا ما مه نیمه ما 10۲ ومبله۷ ۷۱۵۱0 منمعوو 4ص تامنتقجمم. متامهت بط تمعن هچ 7۵0«مطو ف)اناوع۳ فط 1 

فص تقمم 0صه صعمطزمه مط صا مه ومیاله۷ ۱۵09۷۵5۲ مطا 20 رصتقام مط ۵۶ فجعته ممامجعط):۵ه 20 بلحتاصهه مصتعطنامو 
۶ 3690 مضه ,الحگصته: ۵۶ 9۹۵ ,عوهع2۷ ده )92 ٩0۱760‏ فالناوع۴ مطا رووععمعم صمتاهعطانله» فطا تمانه صتقام عطع که ]ماه 
جح بمامم‌لفتامه ۵ جه 02860 0مصصومه ۱۷۵۲۵ فعمتاجبعوهاه ففقط1. .ععقننامد فطا مه مققامععم تم۳/2 111182160 
4 .0۵۳008 ۷۵11020109 مج صمتاه‌طانلق امه 10۲ امممهر فص ۵۶ 1۵۷۵1 قاطاها مه فطا ۵۶ متقصتا ع1صماع2008 
مط ,و۸ ۲۵٩06۵0۷۵۷۰‏ رح 0.34 280 0.35 ۱۷۵۲۵ موه م2)10تاه۷ مه صمتاهتطانل 10۲ متتعاتن ۴۳۷۲9۲ 1دعتادتاماه 
۲ 0.89 280 0.93 220 ,261100 ممتنح‌طنلی عمط 10 0.91 مج 0.94 مج 0متفصتامی مهب متماتی ۱9 لمه ۳ آن داآناوعر 


۶ ۵۱۱626 رعطنهءصنم۴۸ ممتامصصداهعک 4صه ممتاحعتس1 ۵۶ اممصاتهمعن1 ررمووع]۵زظ عاهه0وو۸ 220 )مهن ۲8,۲2۰ -2 20 1 
۰ ,1۲21 ,زه‌ تک رصهتطه 1 ۵۶ الویه۲ نا ردهمعتاموفک لمتنادل 4صه عستانهتتهخ 

(0].26.15 6۵ 02191020 :اتف ممطننه مصت‌طمصومتزمن) -۶) 

صو رده 1 ,برعمامصطهع 1 ۵۶ نمهب نما فده ,همم 0 0۶ ]ممصانجمع1۱ رر0ووع]۳۶۵ -3 
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متحلصته ما عاصاج 7۵ مه هه لاله اوج۲0ه فه۳7 ام0مصه مه عقطا ومتهعن4ص مامتا رمتتهم ممتاهل‌تله۷ معط 
60 ۱۵6 2۲ظ) ٩0‏ رقفطنهه: 01۵20۲۵ م0 0۵۷۵۱0۵۵0 فمامرهع۵ بط مناج صتامتع مظ1. ,و1 تمه« 0صتامع 
م6 26095 اصمصممهصه۱ظ صم2ع1۳1 لصح عویها 0صق1 طا معصقطه فط مصتعتوهمی راهم ۱۷۷/۸ عمط رها 60تمصصتاوه فعباه۷ 
٩30۱۷64 2‏ ومن0ه معط رامعم .تعکننامه ما ققممع اهبع1 ماه 0صنمع ماقاصته زاتهممعج ما ماه م۷۵۲۵ رصتقفام 
تمه فطا ۶ه مصمنعم صماوعب امد 0ص صمط‌نامو . فطع صا مماقم. آوبع1 . مه صتامتع_ ظا وم فتامتامتاومم 
5 0ت)مصصیاوطمع-حع تن ۵۶ ععصعوعيم فطل مه عبتل (رطومطک درم طویعه مه رها ۵۵ رطقطعتیطاه۲ ۵۶ وتمامططم2ع1م) 

۰ 00۵۵۲۵010 0۴ ۲عطاصصباه تمطاوتط عطق 0ص2 م2212 200 ماممه وه طمناه 


«متعصل‌صم) 

۵ ۵061 5۷۷۸۵۸۲ 2060طانلهه مطا حصم 0مصتقاهاه فمله۷ مععقطهع: مطا عقطا 0ع«مطو رنه فنطا گم فلتاوع مط 1" 
6 .165۱118 ۵061 فا ه0۷عمصص دااصهمتصروزه مق هه ماه صنمو ۵ ممتاهانه مومت ت10 قمامصصهتوم آهتمناه 
۶ نامه مها قفا مد وهع۷ )صمع م1 ممصهاهها. ماه0صامع فطا ۵۶ فامممممومم صتقمه مطا ۵۶ فالناوع۲ 
۵ و 10101200 ,مواظ ۷۵۵۲۰ 620 10۲ امممل مو ممتامع1۳0 امه رللحگصته: ۶ ممتافع‌اصا مطا م] عسل معتقطهع۲ 
2 عص1۳0162 رودع۷ مهع07عه فتعامصظ منمانت ممئااتمط 50 م6 30 عنامهاه حصم به۷ وعمتاموع۲ ما7۵ صتا0تع هه 206ناو 
6۵ 22 ۵061 ۲-۱۷۲۵۲۱۳۲,۵۱۷۷-۸۱۷۲۲ ۷۷۵ عطا رلوعصهع 1 .۵0و۲۵ عماج عمط معع/۲عه عمتامعحتعاصا کصه‌نطمزه 
48ص وعمریاموع. مه صامو ‏ مه ممهساه ۵ اممصصمعفصفصه م۲0 101 آ0ما ‏ همهم 2 فه 19860 
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شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی 
5۱۷۸۰۲۷۲۵۵۲۱0۱۷۰ 
(مطالعه موردی: دشت مهاباد) 


امید رجا "- مسعود پارسی نداد" *- مسعود تجریشی ۲ 


تاریخ دریافت: ۱۴۰۰/۱۱/۱۱ 


تاریخ پذیرش: ۱۴۰۰/۱۲/۱۶ 


چکیده 


آب‌های سطحی و زیرزمینی در مقیاس‌های مختلف مکانی یا زمانی باهم در ارتباط می‌باشند. از طرفی منابع آب سطحی و زیرزمینی به‌صورت 
تلفیقی در کشاورزی استفاده می‌شود؛ بنابراین مولفه‌های بیلان آب زیرزمینی و سطحی باید با دقت مناسبی تعیین گردند. در اين مطالعه, تأثیر مقادیر 
تغذیه حاصل از مدل ٩۱۷۸۸۲‏ به‌عنوان یکی از مهم‌ترین مولفه‌های ورودی مدل‌های آب زیرزمینی در شبیه‌سازی تراز سطح ایستابی و مولفه‌های بیلان 
آب زیرزمینی با استفاده از مدل ۱۷۲01(۳۲0۷۷-۱۷۷۲ به‌عنوان یک مدل جامع و یکپارچه منابع آب سطحی و زبرزمینی در دشت مهاباد بررسی شد. 
با استفاده از مقادیر یکنواخت تغذیه و به‌صورت درصدی از بارندگی و کل آب آبیاری در سطح آبخوان بدون توجه به تأثیر تفاوت برداشت و نفوذ عمقی در 
سطح کاربری اراضی در موقعیت‌های مختلف نمی‌توان به‌درستی نوسانات سطح آب زیرزمینی در سطح آبخوان را شبیه‌سازی کرد؛ لذا یکی از مراحل 
مهم تعیین مقادیر آبیاری با توجه به تفاوت کاربری اراضی و تغذیه ناشی از نفوذ عمقی بارندگی و آب ابیاری از منابع آب سطحی و زیرزمینی بود. این 
مهم توسط مدل 5۱۷۸۲ انجام و به‌عنوان شرایط مرزی به مدل ۱۷۲01071,0۷۷-3۷۷۲ تعریف شد. واسنجی و اعتبارسنجی مدل آب زیرزمینی نیز 
به روش‌های سعی و خطا و روش خودکار ۳۳5 انجام شد. دوره شبیه‌سازی به مدت ۱۰ سال از سال آبی ۱۳۸۸-۸۹ تا ۱۳۹۷-۹۸ انجام شد که به 
ترتیب ۶و ۴ سال ۳۹ از ۱۳۸۸-۸۹ تا ۱۳۹۳-۹۴ و ۱۳۹۳-۹۵ تا ۱۳۹۷-۹۸ به‌عنوان دوره واسنجی و اعتبارسنجی در نظر گرفته شدند. با توجه به 
این که الگوی تشکیل مولفه‌های بیلان آب زیرزمینی برای سال‌های آبی مختلف متفاوت است؛ لذا مولفه‌های بیلان آب زیرزمینی برای سه سال معرف 
خشک‌سالی» ترسالی و نرمال نیز بررسی شد. هدایت هیدرولیکی و ضریب ذخیره ابخوان پارامترهای ورودی و واسنجی در مدل ۷۵۲۳۲۷۷۰ 
1 بودند. نتایج واسنجی نشان داد بیش‌ترین میزان هدایت هیدرولیکی و آبدهی ویژه در نواحی جنوبی» مرکزی و شمال شرقی دشت است و در 
حاشیه شمالی و در نزدیکی خروجی دشت به کمترین میزان خود می‌رسد. پس از فرآیند واسنجی, نتایج نشان داد به‌طور متوسط حدود ٩‏ درصد از آب 
پاز نگیو فرمد از مارگ ابا دز بکشن کقاوروی یه ایقوانن تقیقح که ام حاکی آثبرآوره صایت عفن و قایل قیول راد طخ استای سل 
برای هر دو دوره واسنجی و اعتبارسنجی بود. به‌طوری که ننایج معیار 1۷517 برای تراز سطح ایستابی در دوره‌های واسنجی و اعتبارسنجی به ترتیب 
مقدار خطای ۰/۳۵ و ۰1۲۴ متر به دست آمد. همچنین. نتایج معیارهای ۳ و 12 نیز برای دوره واسنجی ۰/۹۴ و ۰/٩۱‏ و برای دوره اعتبارسنجی 
۳ و ۰/۸ برآورد شد که موّید آن است که مدل به‌خوبی واسنجی شده و قادر بوده نوسانات سطح آب زیرزمینی را با دقت مناسبی شبیه‌سازی کند. 
بررسی نتایج اجزای اصلی بیلان آب زیرزمینی برای سه سال آبی ۱۳۹۵-۹۶ (سال خشک)» ۱۳۹۶-۹۷ (سال نرمال) و ۱۳۹۷-۹۸ (سال‌تر) نان داد 
میزان تغذیه ناشی از نفوذ آب بارندگی و جریانات برگشتی آب کشاورزی در هر سه سال متفاوت است. همچنین, فعل‌وانفعالات بین منابع آب سطحی و 
زیرزمینی بین سال‌های آبی از حدود ۳۰ تا ۵۰ میلیون متر مکعب متغیر است که نشان‌دهنده اندرکنش قابل‌توجه بین این منابع است. به‌طور کلی ننایج 


۱ و ۲- به‌ترتیب دانشجوی دکتری و دانشیار گروه مهندسی آبیاری و آبادانی» پردیس کشاورزی و منابع طبیعی» دانشگاه تهران» کرج» ایران 
(#- نویسنده مسئول: 10.26.1۲ ۵) 20 [م51 02 :1ت2ص) 
۳- استاد دانشکده مهندسی عمران, دانشگاه صنعتی شریف تهران» ایران 
۷۷,068 10.22067/15 :101 


۴ شریه آب و خاک. جلد ۰.۳۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


این مطالعه نشان می‌دهد اعمال تغذیه تخمینی مدل ٩۷۷۸‏ همراه با واسنجی هدایت هیدرولیکی و آبدهی ویژه می‌تواند در بهبود برآورد نوسانات 
سطح آب زیرزمینی توسط ۷۷۲(-۷101(۳1,0۱۷ کمک شایانی نماید. نهایتا از مدل تلفیقی می‌توان به‌عنوان یک ابزار کاربردی در تبیین الگوی 
بهره‌برداری مناسب از منابع آب تلفیقی سطحی و زیرزمینی تحت تأثیر سناریوهای مختلف مدیریتی استفاده کرد. 


واژه‌های کلیدی: تراز سطح ایستابی» دشت مهاباه مولفه تغذیه, ۷۷۸۲ ۷0۲۳۲0۷۷-۷۷۲ 


مقدمه 
(2010 .61 61 تا در اغلب مواقع» منابع آب برای کشاورزی 
به‌صورت تلفیقی از منابع آب سطحی و زیرزمینی تامین می‌شود. 
تبادلات تغذیه و برداشت از منابع آب سطحی و زیرزمینی با توجه به 
جانمایی رودخانه‌ها 9 موقعیت چاه‌های بهره‌برداری 9 مداخله انسان در 
مدیریت تأمین آب آبیاری در مناطق مختلف یک چالش منحصربه‌فرد 
در این زمینه است (2019 ,131167 4 6 ۱۷). آب‌های سطحی و 
زیرزمینی در سامانه هیدرولوژی اجزای جدا از هم نیستند و در 
مقیاس‌های مختلف مکانی یا زمانی دارای برهم‌کنشی پویا می‌باشند 
(2016 ,61 ۵ ۱۷۱۵۵۵02 لمج2مک ر.اه ۵ تطمصلمه‌طنهو ر 
217 لذا برای دستیابی به بهره‌برداری بهینه و پایدار از منابع آب» 
مولفه‌های بیلان آب زیرزمینی و سطحی (تغذیه و تخلیه) و پتانسیل 
منابع آب زیرزمینی باید با دقت نسبی برآورد گردد ( .41 6 کصات۸ 
20094 2016 ,.۵1 6۶ 1عطه][۱۷۲۵ :2014 ,اه 6 تتعووهلا ) 
برنامه‌ریزی و مدیریت استفاده پایدار منابع آب به‌ویژه در شرایط 
کم‌ابی» مستلزم شناخت کامل و جامع از وضعیت مبهم و پیچیده‌ی آن 
است. بدیهی است لازمه‌ی این شناخت. در اختبار داشتن ابزاری است 


که بتواند کلیه فرآیندهای در حال انجام در سامانه‌های سطحی و 
با خصوصیات خاص مقیاس زمانی و مکانی که بتواند اثربخشی 
الگوهای بهره‌برداری تلفیقی بر مولفه‌های بیلان آب در محدوده 
مطالعاتی را برآورد کند به‌طورقطع خواهد توانست زیرساخت‌های لازم 
برای حل تناقضات بین عرضه و تقاضای منابع آب را فراهم آورد 
(2018 ,۵ 6۶ حصدل0مطع۱۷۲۵ صحلتهک). 
مدل ۷0۳1/۷۷ به‌عنوان یک مدل استاندارد شبیه‌سازی 

جریان آب زیرزمینی در حالت پایدار و ناپایدار پذیرفته شده است 
(2007 ,۵۵/6۷۰ م2 فصصله2008:۳ رراه 2 ت001ا212272ظ : 
٩606۵۲ ۶ ۵1., 20121‏ :2009 ,۵ ۶۶ 619202). در ین مدل 
مقادیر تغذیه ناشی از بارندگی و جریانات برگشتی آب آبیاری به‌عنوان 
یکی از مهم‌ترین موّلفه‌های ورودی مدل‌های آب‌های زیرزمینی در 
مناطق خشک و نیمه‌خشک به‌عنوان درصدی از بارندگی و کل آب 
آبیاری تعریف می‌شود (:2006 ,.اه ۶ عماصفگ ,عنامع500000 

5 ارزیابی‌ها نشان داده است که عدم قطعیت در تعیین این 
پارامتر, به میزان قابل‌توجهی در نتایج شبیه‌سازی جریان آب زیرزمینی 
تأثیر خواهد داشت. 


صمتاه تج عصعتاه مرح و۲ ااحصنمظ 
تبخیروتعرق بارندگی 
حمنامعنسا ععو۹ 
استفاده از منابع آب 
جهلة تامین مصنا ععمن؟ 


ععققطهع: ۲24( 


ات زیرزمینی ۱ ۱ 


حمت/دعنا ما عصزهه۳ 
پمپاژ از آب زیرزمینی 
جهت تامین آب ابیاری 


رواثاب سطحی 


بت 


آب آبیاری 


1 12221 20 4عاوتامعلا 
۳ جریان جانبی" . 30(:.: ناحیه غیراشباع 
۲ عاه) عاوسلصنهن6 
لته نت سطح آب زیرزمینی 
تبخیر از سطح آب زیرزمینی 
1 5 کر رن 2۶ ممت/هساه5 
۳9 ۳ ععتعطه‌عنل امامت ناحیه اشباع 
تغذیه آبخوان توسط نشت از آبخوان به رودخانه 
رودخانه 


شکل ۱- شمایی از کلیه فرآیندهای سطحی و زیرزمینی 
۹ ۱۱06۵۰0۱0۱ 20 ع6هتای اج 0۶ »تاحصعمداو -1 عساعز۲ 


رجا و همکاران شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی... ۳۵ 


درنتیجه تا زمانی که نرخ دقیق تغذیه به‌عنوان ورودی مدل 
جریان آب زیرزمینی برآورد نسود نمی‌توان با اطمینان رفتار 
طولانی‌مدت یک آبخوان و تأثیرات برداشت از آب‌های زیرزمینی را 
تحت سناریوهای مختلف مدیریتی تخمین زد (2010 ,۷۷۳6۵۵6۲ 

روش‌ها و مدل‌های متنوعی برای تعیین میزان تغذیه آب‌های 
زیرزمینی موجود است. بررسی‌ها نشان داده است که مدل 5۷۷۸۲ 
به‌عنوان یک مدل جامع شبیه‌سازی آب‌های سطحی قادر به برآورد 
تبخیروتعرق واقعی و تغذیه ناشی از بارندگی و آبیاری در مقیاس 
حوضه آبریز و واحدهای پاسخ هیدرولوژیکی در سطح حوضه است 
(:2007 ره ۶ 82جمصه:ز2005۳۴ رفناه‌عآ00ظ502۳ ). از طرفی 
مدل ٩۱۷۸۲‏ در شبیه‌سازی جریانات زیرسطحی دارای محدودیت 
است به‌طوری‌که ماژول آب زیرزمینی آن پارامترهای توزیعی 
هیدروژئولوژی آب‌های زیرزمینی همچون ضریب هدایت هیدرولیکی 
و آبدهی ویژه را در نظر نمی‌گیرد (:2019 ,۵ ۶ نتهزنلظ 61 حصنک 
8 ,.). به‌عبارت‌دیگر فرآیندهای آب زیرزمینی را بیش‌ازحد ساده 
می‌کند و عملکرد ضعیفی دارد. علاوه‌براین» نمی‌توان تأثیر برداشت و 
پمپاژ آب از چاه‌ها را بر شیب هیدرولیکی سطح آب زیرزمینی 
شبیه‌سازی کرد (2010 .۵1 ۵ )2019 .۵ ۵1 صصناط و 

ارزیابی‌ه | نان می‌دهد تلفیق مدل‌ه ای ٩۱۷۸۲‏ و 
۷ می‌توانند با توجه به شرایط هیدرولوژیکی در ابعاد 
مختلف مکانی یک منطقه. امکان شبیه‌سازی دقیق‌تر تبادلات 
جریانات آب‌های سطحی رودخانه‌ها و انهار طبیعی و زیرزمینی را 
فراهم آورد (:2007 .۵1 61 0طمصه:زظ :2005 ,ونامع0ط00٩‏ 
رفط۳۵۳ 240 فامعام‌مطممی.. :2005 مطفتصمت هه صنا 
0 با استفاده هم‌زمان از دو مدل برآورد نفوذ حاصل از رواناب 
سطحی و جریان جانبی در خاک به شبکه‌ی جریان و امکان بررسی 
بیلان آب زیرزمینی و برآورد نوسانات سطح ایستابی فراهم می‌گردد 
(:1 1 , .۸1 6۶ ۷۵۲ 200 ۱0۳/۱22020 :2016 ,.آه 6۶ ۳116۷ 

09 #7 حصناطن ‏ :2016 ,۵01۳00۳0120 2). همچنین. در 

صورت استفاده هم‌زمان نه‌تنها محدودیت‌های دو مدل (برآورد تغذیه 
در مدل ۷100۳۲0۷ و جریانات زیرسطحی در 57۷۵۸۲) را می‌توان 
پوئش داد بلکه خصوصیات زمانی و مکانی منطقه موردنظر را به‌طور 
مناسبی منعکس کرد (2008 و 6۶ لک ررماتوظ مه عانهظ ۲‏ 
75 2018 ,.ه 6 :۷۷۵ ). به‌عبارت‌دیگر تغذیه حاصل از 
مدیریت‌های مختلف آبیاری شامل تلفات نفوذ عمقی ناشی از تأمین 
مقادیر متفاوت آب آبیاری و مقادیر تغذیه ناشی از نفوذ بارش و نشت 
از رودخانه‌ها در سال‌های مختلف آبی و در پهنه‌های مختلف دشت 


۱۳ ۱ 


حاصل از مدل 81۷۸۲ به‌عنوان ورودی مدل آب زیرزمینی 
۷ مورداستفاده قرار می‌گیرد. مزیت مهم دیگر مدل 
تلفیقی» بررسی موّلفه‌های بیلان و اندرکنش آب‌های سطحی و 
زیرزمینی در مقیاس‌های زمانی و مکانی مولفه‌های بیلان است. 
درعین‌حال می‌توان تغییرات سطح ایستابی و شیب هیدرولیکی آب 
زیرزمینی تحت الگوهای مختلف برداشت و مدیریت بهرهبرداری» 
میزان نرخ برداشت از چاه‌های بهره‌برداری از آبخوان تحت سناریوهای 
مختلف را مدیریت کرد (:2012 ,.۵1 ۶۶ 6061۵۲ ر.اه ۶ 17207 
90 ,2556۲ :2012). 

از یک‌سو سیر گستره نوع فعالیت‌های بهره‌برداری از منابع آب و 
عدم توجه کافی به برنامه‌ریزی و مدیریت تلفیقی منابع آب و از طرف 
دیگر عدم توازن مکانی و زمانی توزیع منابع آب و محدودیت‌های 
استفاده از آب‌های سطحی, ضرورت تحقیق در خصوص بهره‌برداری 
تلفیقی از منابع آب سطحی و زیرزمینی را ایجاب می‌نماید 
(:20186 .اه ۶ تعطهاز]20191۷ ,مات 220 ۷۷۵ ). با توجه به 
گسترش مدل‌های عددی طی دو دهه گذشته, هنوز هم کمبود 
استفاده از مدل‌های جامع برای ارزیابی و مدیریت همه‌جانبه و 


یکپارچه منابع آب سطحی و زیرزمینی وجود دارد ( 127020 
8 ( ۵ ۵ ۱۷۱۵۵2 رراه ۵۶ صهل0هطوه۷]۵ 2۵8ع ٍ 


6 به‌ویژه اتصال مدل‌های آب سطحی و زیرزمینی به‌طوری که 
زیرزمینی» بلکه تأثیر سامانه آب زیرزمینی در منطقه غیراشباع و 
به‌ویژه منطقه توسعه ريشه را نیز مدنظر قرار دهد ( ,4 ۶ ۳۲116۷ 
20164 6۶ ۷۷۵ :2019 ,۲2116۷ 200 ۷۷۵ :2019 ,.۵1 ۶ صصتال 

58 ,.0). از طرفی در دهه‌های اخیر تأمین آب موردنیاز از منابع 
زیرزمینی منجر به توسعه استراتژی‌های مدیریتی توآمان بهره‌برداری 
از منابع اب سطحی و زیرزمینی شده است. به‌طور مشخص در حوضه 
دریاچه ارومیه که به‌شدت تحت تأثیر بهره‌برداری غیرمتعارف از منابع 
آب قرارگرفته. این مدل‌ها می‌تواند برای حداکتر سازی تأمین حقابه 
دریاچه ارومیه با توجه به الگوی تأمین نیاز آبیاری اراضی کشاورزی از 
منابع تلفیقی مورداستفاده قرار گیرد. این مطالعه تلاش دارد به بررسی 
تأثیر مقادیر تغذیه ناشی از نفوذ عمقی حاصل از تأمین توآمان آب 
آبیاری از منابع آب سطحی و زیرزمینی و نیز نفوذ ناشی از بارندگی 
حاصل از مدل ٩۱۷۸۲‏ در شبیه‌سازی آب زیرزمینی با استفاده از مدل 
۱1۳0110۷۷۲ به‌عنوان یکی از پارامترهای بسیار مهم و 
تأثیرگذار در تعیین مولفه‌های بیلان آب زیرزمینی در دشت مهاباد 
بپردازد. مولفه‌های بیلان آب زیرزمینی در سال‌های آبی متفاوت 
به‌شدت تحت تآثیر تغذیه و تخلیه قرار می‌گیرد لذا موّلفه‌های بیلان 
آب زیرزمینی در منطقه موردمطالعه در سه سال آبی خشک. نرمال و 
مرطوب نیز تجزیه‌وتحلیل خواهد شد. 


۶ شریه آب و خاک. جلد ۳۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


ایستگاه هیدر ومتری گرد یعقوب 00اه با 
سد مهایاد 
سد انحرافی بوسف گندی 1۵17 ۱0 


۲ 0 9 ۹ 4 


0۶ 


ت_ 3 0[ 


ماد تست متداییه ااب ست. ام سس 
سطح آب زیرزمینی سطح زمین سنگ پستر 
۳ مرو 
مود ۳0 
3 
مود 
2 ووعه 
مود « ِ 
وود سد ز 1 
مهد 33 
1 محد 3 
۳ 
3 مود 
نك 
265 فت 
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(۳) مها 
فاصله (متر) 


0 


شکل ۲- موقعیت شبکه آبراهه. سد مهاباد. سد انحرافی یوسف کندی, ایستگاه هیدرومتری گرد بعقوب. آبخوان مهاباد در حوضه آبریز ارومیه و 


وضعیت تراز سطح ایستابی در آبخوان مهاباد در طول دوره آماری ۱۳۷۵-۹۸ 
و6)24101 0۳06016 و حابمامهو0 نع رصق ممتممه لومعم همان ۷ مه 0هصامصما باد۵ ۵0۲ عصصمع‌تمای ۵ ممتای م1۵ -2 ۲1۵0۲۵ 
1996-09 متس گنود ۱۷۸۵۵0۵20 مر (۱6۲۷۵ ۲۵6۵ 200 مصتعوظ متصصرنا صز وه ۷۱22020 


مواد و روش‌ها 
منطقه موردمطالعه 

شوه لیا ره بیان او سای ما سم ۱۵۶ 
کیلومترمربع» هشتمین زیرحوضه از حوضه آبریز درياچه ارومیه است 
که در شمال غربی ایران و جنوب غربی حوضه آبریز دریاچه ارومیه 
واقع‌شده است. موقعیت این محدوده در کشور و حوضه آبریز دریاچه 
ارومیه در شکل ۲ نشان داده‌شده است. در این محدوده می‌توان به 
سد مپاباد اشاره کرد که از سازه‌های مهم در منطقه است و شبکه 
آبیاری و زهکشی مهاباد در پایین‌دست آن قرار دارد. در این محدوده» 
وسعت آبخوان مهاباد حدود ۲۵۳ کیلومترمربع است. به لحاظ اقلیمی 
مطابق با روش‌های دمارتن و آمبرژه دارای آب‌وهوای نیمه‌خشک و 
متونط پارندگی‌یر نانی آمازباندگی ایسشگاه سیتوشیک: با در 
طول دوره آماری ۱۳۶۸-۱۳۶۹ تا ۱۳۹۷-۱۳۹۸ سالانه ۴۰۲ میلی‌متر 
کرازان شفه سیف یل شوه اما دکزطی مش رب یکی اه 
از تشت بر اساس آمار و اطلاعات جمع‌آوری‌شده ایستگاه سینوپتیک 
مهاباد برابر یا ۱۸۶۰ میلی‌متر است. متوسط حداکفر و حداقل ماهائه 
دما با میزان ۳۳/۷ و ۲/۲- درجه سانتی‌گراد به ترتیب در ماه‌های 
مرداد و بهمن اتفاق افتاده است. روند تغییرات سطح آب زیرزمینی در 
طول ۲۳ سال (۱۳۷۵-۹۸) نشان می‌دهد سطح آب زیرزمینی به 
میزان ۱/۵ متر افت داشته است که متوسط سالانه آن برابر با حدود ۶ 
سانتی‌متر در سال است. افت سطح آب زیرزمینی عمدتاً به دلیل 
برداشت آب برای تأمین آب آبیاری و تبخیر از آب زیرزمینی به دلیل 
بالا بودن سطح آب زیرزمینی در مناطق شمالی و مرکزی دشت است 
(شکل ۲). بررسی‌ها براساس مقطع طولی از سطح آبخوان و با توجه 
به تراز سطح زمین» آب زیرزمینی و سنگ‌بستر نشان می‌دهد ضخامت 
سنگ‌بستر از جنوب‌غربی به شمال شرقی بین ۲۱ تا ۵۳ متر و تراز آب 


زیرزمینی از ۱/۵ متر در نواحی شمالی تا ۶/۵ متر در نواحی جنوبی 
آبخوان متغیر است. به‌طور کلی روند اصلی تراز سطح آب دشت در 
سال‌های اخیر متعادل بوده و افت سالانه مشاهده‌شده در مقایسه با 
بسیاری از دشت‌های کشور ناچیز است. نوسانات سطح آب در دوره 
بیلان از یک منحنی سینوسی پیروی کرده و به‌صورت شاخه‌های 
صعودی و نزولی تقسیم می‌گردد. از شروع سال آبی در مهر تا 
اردیبهشت ماه شاخه صعودی بوده که به عنوان دوره مرطوب قلمداد 
می‌گردد. از اردیبهشت تا شهریور ماه منحنی نزولی ادامه می‌یابد که 
به‌عنوان دوره خشک هیدروگراف قلمداد می‌گردد. به همین دلیل 
هیدروگراف ارائه شده دارای فراز و فرودهای فصلی ناشی از تغذیه در 
طول دوره مرطوب و پمپاژ از آب زیرزمینی در فصل‌های گرم و عدم 
تغذیه آن از طریق بارش در دوره خشک می‌باشد. 


تئوری مدل ۷۲0۲۲۲۸۵۷۲۷ 

مدل ۷]01(771.01#7 اولین بار در سازمان زمین‌شناسی آمریکا! 
(۲505]) در سال ۱۹۸۲ تهیه و از آن زمان تاکنون نسخه‌های فراوانی 
از آن منتشرشده است (1995 ,.۸1 ۰۶ ۷101>66). این کد نرم‌افزاری با 
توجه به قابلیت‌های بالا و داشتن پکیج‌های مختلف که قابلیت 
شبیه‌سازی اجزاء مختلف تأثیرگذار بر سامانه آبی را دارنده یک مدل 
کامل در مدل‌سازی عددی یک سامانه اشباع و تشکیلات آبرفتی 
است. مدل ۷۲۵۲۱۳۲۸۵۷۷[ یک مدل سه‌بعدای جریان اشباع» تفاضل 
محدود و 131001-00016760 است که جریان آب زیرزمینی در محیط 
متخلخل را به‌صورت سه‌بعدی شبیه‌سازی می‌کند ( .اه 2 1۷0۲66 
5 برنامه‌ی مذکور به دلیل قابلیت و اعتبارسنجی گسترده 
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به‌عنوان یک مدل استاندارد جریان آب زیرزمینی در سراسر جهان 
پذیرفته شده و در شبیه‌سازی و تجزیه‌وتحلیل‌های حالت پایدار و 
ناپایدار جریان آب در سامانه‌های مختلف استفاده می‌شود ( 42صهطن5 
20094 ,.۵1 6 ,686۷ظ 0ص حصصله۲ :2008 .اه 61 011ا2۲2272ظ۳ 
2 ره ۶ 56062۲ :2007). مدل ۷01۲1/۵۷۷ مدلی 


فیزیک محور و کاملاً توزیمی است که معادله دارسی را با بیلان جرم 
برای جریان زیرسطحی ترکیب و از معادله حاکم بر جریان آب 
زیرزمینی برای شبیه‌سازی جریان استفاده می‌کند ( .اه 6 10 
8 معادله‌ی دیفرانسیل جزئی ذیل که از ترکیب شکل سه‌بعدی 
معادله دارسی و معادله پیوستگی حاصل می‌شود. معادله حاکم بر 
جریان آب زیرزمینی در محیط اشباع است: 


ث 5 )+ ِ )+ 5 ) 
)۱ 02 22 2 رد۵ 2۲۲ 0 رن ۳ 02 


که در آن »کل ورگ و یک به ترتیب مقادیر هدایت هیدرولیکی در 
جهت‌های »* لو 2 ,5 ضریب ذخیره ویژه (41/1 7 بار آبی» /زمان 
است؛ و ۲۷ فلاکس جریان حجمی در واحد حجم (/1) است که 
برای تخلیه منفی و برای تغذیه مثبت است. اين معادله با توجه به 
شرایط آبخوان موردمطالعه مبنی بر همگن یا غیر همگن بودن 
همروند یا غیرهمروند بودن و ماندگار یا غیر ماندگار بودن قابل 
ساده‌سازی است. 

این مدل برای شبیه‌سازی جریان آب زیرزمینی شامل برنامه 
اصلی و تعداد زیادی برنامه‌های فرعی بوده که این برنامه‌های فرعی 
در تعدادی نرم‌افزار مستقل دسته‌بندی شده‌اند. هر بسته نرم‌افزار برای 
شبیه‌سازی یکی از سامانه‌های هیدرولوژیکی نظیر تغذیه از رودخانه» 
فیه #وسا اهقدیدافت آب تسط اه برانضل سادلات 
خطی جریان با یک روش خاص به کار می‌رود. به‌طوریکه هرکدام از 
پی‌های محاسای و هیترواوژیگی به‌گونهای ایضاه شته‌اند که 
مستقل از یکدیگر عمل می‌کنند و به دلیل وجود این ساختار بخش‌ها و 
گزینه‌های جدیدی را می‌توان به مدل افزود بدون آنکه نیاز به تغییر در 
بخش‌های موجود باشد. بسته رودخانه" با ساده‌سازی نحوه‌ی ارتباط 
بین آبخوان و رودخانه و میزان جریان مبادله شده بین آن دو را 
معطوف به دو عامل اختلاف هد و ضریب هدایت بستر رودخانه 
می‌نماید (2016 ,اه ۵۶ 92220ا)). بسته‌های شبیه‌ساز 
تبخیروتعرق در ,۷10۵1۳۲ از روش تقریب خطی برای محاسبه 
میزان تبخیروتعرق استفاده می‌کنند. در این روش نرخ تبخیروتعرق 
روی سطح زمین و حداکثر ارتفاعی که در آن تبخیروتعرق از آب 
زیرزمینی رخ می‌دهد توسط کاربر مشخص می‌گردد و بسته شبیه‌ساز 
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تبخیروتعرق با استفاده از داده‌های دریافتی از کاربر یک تقریب خطی 
ساده و يا یک تقریب دو یا چندخطی از پروفیل نرخ تبخیروتعرق در 
عمق زمین ایجاد نموده و میزان آب خارج شده از آبخوان را محاسبه 
می‌نماید. 

یکی از نسخه‌های جدید مدل ۷۲0۲3۳],۷۷ نسخه 
2-۲ ۷۵۳۳1,0۱۷ است که در اين تحقیق از آن برای ساخت 
مدل آب زیرزمینی دشت مهاباد استفاده شد. در مدل -۱10۳۲,0۷۷ 
1 از فرمول نیوتن برای حل مشکلات مربوط به خشک شدن 
سلول‌ها و دوباره خیس شدن غیرخطی در معادله زیرزمینی جریان 
محصور نشده استفاده می‌کند. این قابلیت در نسخه‌های قبلی مدل 
۷0۷ وجود نداشت و جز نقاط ضعف آن‌ها بود. این صدل 
باید با بسته (۳۷۷) 0012ع0506270-۷۷/] استفاده شود. این 
بسته برای محاسبه قابلیت هدایت داخل سلولی از روشی متفاوت از 
روش (0۳ظ) ۳1۵۳ 31001-0670060760 استفاده می‌کند. به‌طور 
خلاصه جنبه‌های عملی مهم مدل ۱۵۳0۳10۱۷۷۷۲ که 
متفاوت از مدل ۷6010۳1017-2005 است شامل موارد زیر است: 
۱) از عطناطعنه۷۷-صدع‌تاوونا برای محاسبه هدایت افقی آبخوان- 
های غیر محصور استفاده می‌کند» ۲) تمام سلول‌های بار متغیر فعال 
در شروع شبیه‌سازی» در سرتاسر شبیه‌سازی فعال باقی می‌مانند» ۳) 
هدایت افقی برای شرایط غیر محصور و ناپایدار در طول سلول‌های تر 
و خشک تصحیح‌شده است» ۴) فرمول‌بندی ذخیره برای تغییرات 
ذخیره در طول سلول‌های خشک و تر و در طی شرایط بین محصور و 
غیر محصور تصحیح‌شده است و ۵) فرمول‌بندی ذخیره برای شرایط 
غیر محصور تصحیح‌شده است به‌طوری که تغیبری در ذخیره به علت 
تغییرات بار در زیر کف سلول اتفاق نمی‌افتد ( .۵1 6 ۲15۷0۲6۲ 
2011 


مدل مفهومی منطقه مورد مطالعه 

تهیه مدل مفهومی ساده‌ترین راه برای شناسایی و تبیین وضعیت 
ابخوان است. یک روش گام‌به‌گام مرحله‌ای را برای تهیه‌ی مدل 
مفهومی ۷0](۳۲,۵۷۷-۷۷/۲ در شکل ۲ نشان داده شده است. 
برای تهیه مدل کمی و شبیه‌سازی جربان آبخوان مهاباد در مدل 
۲ 0۷۷ ,۷101۳1 نیاز به ورود یک سری اطلاعات هواشناسی؛ 
هیدرولوژی» ژئوفیزیک» زمین‌شناسی و هیدروژئولوژی از منطقه است. 
به‌طور مشسخص در مدل‌سازی آب زیرزمینی با استفاده از مدل 
۷10۳10۷۷۲ در نرم‌افزار 0115 پوشش‌های " مختلفی 
شامل رقوم سطح آب در چاه‌های مشاهده‌ای» دبی چاه‌های 
بهره‌برداری» اطلاعات هیدرومتری» نوع مرزهای محدوده ابخوان» 
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ورودی و خروجی‌های آب‌های سطحی و زیرزمینی» سنگ کف لایه 
شناسی آبخوان. خصوصیات هیدرولوژیکی و هیدرودینامیکی ابخوان 
ویژه) و جانمایی چاه‌های بهره‌برداری و مشاهده‌ای برای هر سلول 
وارد شد. 


داده‌های مورداستفاده در مدل 


یکی از مهم‌ترین پارامتر ورودی نقشه ۳۷ مناسب است. 
آگاهی از دقت نقشه‌های توپوگرافی در مدل‌سازی آب زیرزمینی دارای 
اهمیت بسیار زبادی است. در اين مطالعه از نقشه 91*11۷ در مقیاس 
۰ متر در تهیه تمامی نقشه‌های آب زیرزمینی و هندسه آبخوان 
مهاباد مانند مرز آبخوان, تراز آب زیرزمینی» سنگ کف و غیره استفاده 
شد (شکل ۴۵). به‌طورکلی» شیب عمومی دشت از جنوب غربی به 
شمال شرقی و متوسط شیب دشت تقریباًبرابر هفت در هزار است. با 
توجه داده‌های در دسترسء بر اساس نقشه‌های زمین‌شناسی و 
توپوگرافی منطقه, نقشه هم‌تراز آب زیرزمینی» موقعیت گسل‌ها در 
منطقه و نیز وضعیت پراکندگی چاه‌های حفرشده در منطقه» آبخوان 


مهاباد به‌صورت تک لایه در نظر گرفته شد. چاه‌های پیزومتری نیز به 
مدل معرفی شد. شبکه چاه‌های پیزومتری به‌منظور اندازه‌گیری سطح 
آب زیرزمینی آبخوان آبرفتی از سال ۱۳۶۹ در دشت مهاباد آغازشده 
است. تعداد چاه‌های مشاهده‌ای دشت مپاباد در سال ۰۱۳۸۰ ۲۰ حلقه 
بوده و در سال‌های بعد به‌تدریج بر تعداد آن‌ها به‌منظور تکمیل شبکه 
چاه‌های مشاهده‌ای افزوده‌ شده است. به‌طوری که در سال ۱۲۹۰ تعداد 
۲ حلقه و درنهایت در سال ۱۳۹۷ نیز تعداد ۵ حلقه چاه حفره شده 
است. طیق بررسی‌های انجام‌شده» از آمار و اطلاعات ۲۲ حلقه چاه 
پیزومتر که سطح آب زیرزمینی به‌صورت ماهانه در آن‌ها اندازه‌گیری 
می‌شود و دارای شرایط آماری مناسب در طول دوره شبیه‌سازی بود در 
مدل‌سازی استفاده شد (شکل ۴0). لازم به ذکر است داده‌های ارتفاع 
آب در پیزومترها در فرآیند مدل‌سازی تأثیر ندارنده اما به‌عنوان داده‌ی 
مشاهده‌ای در ارزیابی نتایج و واسنجی مدل از آن‌ها استفاده می‌شود. 

برای شبیه‌سازی دقیق اثر رودخانه‌ها در مدل. اطلاعات مربوط به 
محل ایستگاه هیدرومتری واقع در دشت (ایستگاه گردیعقوب) از 
سازمان آب منطقه‌ای آذربایجان غربی اخذ گردید. 


چاه‌های بهره‌برداری و مشاهده‌ای 


۱ 
له 


اطلاعات هیدرومتری و مدل رقومی 
ارتفاع 


0 ها م0 »0096۱ 1 
1 س«مذ) «ع۱ع 1و1 


اطلاعات خصوصیات هید رودینامیک 
آ 


اس هرهب وا کنو 
هام یماهس 


اطلاعات مربوط به زهکش و رودخانه 


۱ 


مقادیر تغذیه حاصل از مدل 


| 
من 5۱۲۵۲ ع۱ 


جمع‌آوری داده‌های ورودی به مدل 
0/0( 


۵۶ 0 هام )دورد وهای 
امعم ۱۱۵۵/۵۱۰۱ 


تنظیمات شرایط مرزی گام زماتی و اجرای مدل 


سا 2۵0 ,موی صن) بکع نامک صمتانلصم 130070021 


۱ ۵4 


خروجی مدل 


1۲۶ ۱۵6۱ ۸ 


سطح ایستابی و بیلان آب زیرزمیتی 


ماو وی ۵۳۵0 ۱2۵۱6 ۱۷۵/۵۰ 


واسنجی و اعتیارسنجی مدل براساس پارامتر سطح ایستابی با استفاده از آزمون سعی و خطا و نرم‌افزار ۳۳5۲ 
او عصنهه 09946 ماه ماج 0عقون ۵06 ۵6 گم سمتوتاه۲ 990 ممتاورطناو 
۳۴۹۲ ره بمنا نا ۵0 2۳0 


شکل ۳- داده‌های مورد استفاده در ساختار مدل‌سازی ۷۲۷-۱۲۷۲۷۲ ,5۷۷۸۲-۷۲01۲ 
عصنام0 ۳۵ ۱۲۷-۱۲۲۲ ,]5۱۲۷۲۸۰۲-۱۵۲ رز 60ج 1۵6۵ -3 عع۲1 
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البته به دلیل کم بودن تعداد ایستگاه‌های هیدرومتری که این 
اندازه‌گیری‌ها در ایستگاه مذکور انجام‌شده از دفت در مدل‌سازی 
کاسته می‌شود. برای تقویت اطلاعات ورودی از ساير اطلاعات 
همچون اطلاعات خروجی برآورد شده از مدل ٩۷۷۸۲‏ شامل 
خصوصیات شیب پهنا و عمق رودخانه اصلی مهاباد و کلیه آبراهه‌ها و 
اطلاعات گزارش‌شده در مقالات علمی (2015 ,۱2۷61) استفاده 
شد. شبکه رودخانه مورد استفاده در فرآیند مدل‌سازی در شکل ۴0 
نشان داده‌شده است. 

جریان آب زیرزمینی ورودی یکی از مهم‌ترین عوامل مهم تغذیه 
آبخوان است. لذا جبهه‌های ورودی و خروجی آبخوان مهاباد به‌عنوان 
شرایط مرزی با استفاده از نقشه تراز سطح ایستابی» خطوط هم‌تراز آب 
زیرزمینی و جهت جریان آب زیرزمینی مشخص شد. پس از بررسی 
خطوط هم‌تراز آب زیرزمینی مشخص شد آبخوان مهاباد دارای سه 
مقطع ورودی» یک مقطع خروجی و دو جبهه نفوذناپذیر است (شکل 
0 یکی از مهم‌ترین منبع آب موجود در دشت مهاباده چاه‌های 
بهره‌برداری جهت مصارف شرب صنعت» کشاورزی و سایر مصارف 
است. میزان برداشت از آب‌های زیرزمینی در سال‌های ۱۳۷۵ ۱۳۸۷ 
۰ و ۱۳۹۷ آماربرداری شده است. بررسی‌ها نشان داد به‌طور 
مشخص در طی دو دهه اخیر به‌تدریج تعداد چاه‌های حفرشده 
افزایش یافته و میزان بهره‌برداری از چاه‌ها در کل دوره آماری افزایش 
چشمگیر داشته است. با توجه به بانک اطلاعاتی آمار مطالعات پایه 
منابع آب» میزان برداشت از آب زیرزمینی در سطح دشت از ۱۲/۹۰ 
میلیون مترمکعب در سال آبی ۱۳۸۹-۹۰ به ۵۵/۲۰ میلیون مترمکعب 
در سال ۱۳۹۷-۹۸ (بیش از چهار برابر) از طریق افزایش تعداد 
چاه‌های پمپاژ از ۲۳۱۲ به ۳۴۲۶ افزایش‌یافته است. موقعیت چاه‌های 
بهره‌برداری در آبخوان مهاباد در شکل ۳4 نشان داده‌شده است که 
داده‌های برداشت از آن‌ها در طی فرآیند مدل‌سازی به مدل تعریف 

جهت برآورد رقوم سنگ کف اطلاعات برداشتی مطالعات 


ژئوفیزیک» نقشه‌های زمین‌شناسی, لوگ چاه‌های پیزومتری و 
داده‌های لوگ حفاری چاه‌های اکتشافی و بهره‌برداری از سازمان‌های 
مربوطه موردنیاز است. پس از بررسی‌های لازم از اطلاعات ٩۳‏ حلقه 
چاه جهت رسم رقوم سنگ‌بستر نسبت به تراز سطح زمین در نرم‌افزار 
5 ۸ استفاده شد (شکل ۴۵). به‌طور کلی رقوم سنگ کف آبخوان 
از ۱۲۳۲ تا ۱۲۸۴ متر متغیر است. آبرفت سیلاب‌دشت شامل رسوبات 
آپیسیتا نفوذپ ذیر از شن و رس است ( عصتادعهمت ۷۷۵16 
4 راصمطع6۱0 ۲۱۱۵۷ ماطفمته)منی 24 عتععصتع۳). ضخامت 
سنگ پستر» متشکل از سنگ‌آهک کرتاسه و حدود ۱۶-۶۰ متر زیر 
موم فا شرف 

ضرایب هیدرودینامیک مهم‌ترین عوامل کنترل‌کننده حرکت 
سیال در محیط متخلخل می‌باشند. ضرایب هدایت هیدرولیکی (6 
قابلیت انتقال (7) و آبدهی ویژه (ر5) از دسته پاامترهایی می‌باشند 
که به میزان سیال و تغییرات آن مرتبط است. در اين مطالعه پس از 
بررسی‌های لازم از منابع اطلاعاتی متفاوت از قبیل آزمايش پمپاژ در 
چاه‌های اکتشافی» پیزومتری» بهره‌برداری و بر اساس دانه‌بندی نمونه 
خاک (1996 ,۷۷1650 16 20 1(2۷15) برداشت‌شده در چاه‌های 
بهره‌برداری برای تعیین ضرایب هیدرودینامیک استفاده شد. پهنه‌بندی 
نقشه قابلیت انتقال نشان داد با توجه به ضخامت کم آبرفت و 
کوهستانی بودن منطقه ضریب قابلیت انتقال این دشت دارای ارقام 
پائینی بوده و از ۲۵۰ تا ٩۳۰‏ مترمربع در روز متفیر است. به‌طوری که با 
عبور اراضی جنوب و جنوب غربی به سمت مرکز دشت میزان آن 
کاهش يافته و در مرکز دشت به کمتر از ۵۰۰ مترمربع در روز می‌رسد. 
در نواحی شمالی دشت میزان قابلیت انتقال بین ۲۵۰ تا ۴۰۰ متر مربع 
در روز متغیر است. نتایج مطالعات مهندسین مشاور آب و توسعه پایدار 
| عاطمصتفاعیی مجح قععصتعص . مصلنمجی . هاه۱۷۷ 
4 ,۵۱00۵۳:60 ۲(۵۷) و ولیزادگان و یزدانیناه ( ,بطهعع۷21:220 
8 ,فجوم02 ۷۵2 200) نبز نشان داده است میزان قابلیت 


انتقال بین ۲۵۰ تا ۹۵۰ مترمربع در روز متغیر است. 
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کانال آبپاری +0 ۱۳۳00 چچچچمن وه 


45 45 45 45 45 45 
شکل - ه) مدل رقومی ارتفاع (12۷1) با مقیاس ۲۰ متر 0) موقعیت چاه‌های پیزومتری و شبکه رودخانه» ») مرزهای ورودی خروجی. نفوذناپذیر 
و خطوط هم‌تراز آب زیرزمینی» 0) موقعیت چاه‌هاء ع) مقادیر سنگ‌بستر آبخوان» ؟) مقادیر هدایت هیدرولیکی» ب) موقعیت کانال آبیاری» سد و نقاط 
تحویل آب در سطح شبکه و ) نقشه کاربری اراضی 
00166 و116 0 وواله ۲ 0۱620۳060۵ 01 «منای‌مصا (0 رعماعصه 30 ۵۶ مهو وه ۷۵۵ (۱۸۱۷) [06مصط معا اهانعز۱ (و -4 و۲ 


تاه۲۱۲0 ( مه 06000 نت۸ عماج ۷ ره وعااه ۳ 0۲ «متاج)منبا 0 وعه‌صنا موی 0۲هام 2۳00 ومم021صنوها ماصامعصصه مر 
۵ 0۹۵ ۱۵۱0 (ظ 20 ۱۵۲۱۲۷۵0۲ 6ظ) اج فصمهر متام تماه ۷ 2۳00 صهل ماقصی ممتامفت اه جمتای‌منا رم رممماو تام صمصمم 


رجا و همکاران شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی... ۴۱ 


در اين مطالعه از اطلاعات آزمایش پمپاژ چاه‌های اکتشافی و 
بهره‌برداری برای تعیین هدایت هیدرولیکی و آبدهی ویژه استفاده شد. 
داده‌های حاصل از آزمایش پمپاژ نشان داد میزان هدایت هیدرولیکی 
از ۲/۵ تا ۱۲ متر در روز متغیر است (شکل ۳۶). به‌طورکلی نواحی 
شمالی و شمال غربی آبخوان کم‌ترین مقدار و نواحی جنوبی و جنوب 
غربی بیش‌ترین مقدار را نشان می‌دهند. نزدیک خروجی ابخوان 
میزان هدایت هیدرولیکی به دلیل وجود رسوبات دانه‌ریز کاهش 
می‌یابد. همچنین» میزان آبدهی ویژه در سطح منطقه از ۰/۴ تا ۳/۶ 
مترمربع در روز متغیر است. بررسی مطالعات پیشین به‌منظور تعیین 
آبدهی ویژه نیز نشان می‌دهد مقدار آن بین ۰۱۳ الی ۱ درصد است 
| عاطمصتفاونی 4حصه قععصتعص . عصتانمصم . ماه۱۷۷ 
4 ,1(۵۷۵۱۵۵۳۵). 

محدوده تغییرات ضریب هدایت کف رودخانه‌ها برای سه بازه 
ورودی به دشت. میانی و پایاب برای رودخانه مهاباد بين ۲ تا ۱۵ متر 
در روز متغیر بود که در مرحله واسنجی مدل جهت تخمین اولیه و 
درنهایت تعیین مقدار دقیق ضریب هدایت کف رودخانه استفاده شد 
(2015 رتتهرردلا مصده فععصنعمظ ممتانممم ماه ۳۷ 

4 راطعصصمم1(۵۲۵1 ماطفصته)کنام) درک فرآیندهای تغذیه آب 
زیرزمینی برای توسعه منابع آب در مناطق نیمه‌خشک و خشک مهم 
است. در این مطالعه» در مدل‌سازی آب زیرزمینی از مقادیر تغذیه 
حاصل از مدل ٩۷۷۸۲‏ استفاده شد (2021 ,.۸1 6 6212 میزان 
تغذیه (نفوذ عمقی) ناشی از تلفات آب آبیاری و نفوذ از بارندگی از 
پروفیل غیراشباع خاک حاصل می‌شود. مدل 519۸ فرض می‌کند 
که آب آبیاری از بالاترین لایه به پایین‌ترین لایه خاک تا ظرفیت 
زراعی پر می‌شود و هنگامی‌که ذخیره رطوبتی خاک در یکلابه بیش 
از ظرفیت زراعی باشد» نفوذ عمقی رخ می‌دهد. در پرونده‌های خروجی 
مدل در قسمت 0۱00067 پارامتر 0۷۷_101160) که بیانگر تغذیه 
به آبخوان است. استخراج و به مدل معرفی شد. به‌طورکلی میزان 
تغذیه سالانه واردشده به مدل از ۱۷ تا ۲۶۲ میلی‌متر متغیر است. 
اگرچه به‌طور قاطع نمی‌توان الگوی مکانی مشخصی را برای تغذیه در 
سطح همه کاربری اراضی و در کل سطح دشت اذعان داشت. اما 
نتایج شبیه‌سازی حاصل از مدل ٩۷۷۸۲‏ نشان داد که مقدار 
قابل‌توجهی تغذیه در سطح اراضی زراعی و باغی و در مرکز و خروجی 
دشت (واقع در شمال غربی) و نیز در مناطق واقع در امتداد کانال 
اصلی رودخانه وجود دارد. موقعیت کانال‌های آبیاری به همراه نقاط 
کاربری اراضی استفاده شده در مدل ۶۷۷۸۲ در شکل ۴ نشان داده 
شده است. اطلاعات تکمیلی در خصوص نتایج مدل ٩۷۷۸۲‏ در 
مطالعه رجا و همکاران (2021 ,.۵1 ۵ 2[2) منتشر شده است. 
خشکسالی پدیده‌ای پیچیده است که می‌توان آن ر خستط با 


شدت. مدت زمان و میزان آن مشخص کرد. از بین این سه بعده 
شدت خشکسالی عامل اصلی می‌باشد که می‌تواند برای تجزیه و 
تحلیل خشکسالی استفاده شود. استانداردها برای بیان خشکسالی 
مبتنی بر اقلیم‌شناسی منطقه است. به طور معمول, اندازه‌گیری‌های 
هواشناسی اولین شاخص‌های خشکسالی می‌باشند. خشکسالی 
هواشناسی را می‌توان در نتیجه کاهش و یا عدم بارش مشاهده کرد و 
آن را در یک منطقه برای یک دوره زمانی تعریف کرد. ین شرایط 
معمولاً برای هر منطقه متفاوت است؛ زیرا شرایط جوی در مناطق 
مختلف از نظر مکان و زمان بسیار متفیر می‌باشد. خشکسالی 
هواشناسی با توجه به شاخص‌های تعریف آن بر اساس داده‌های آب و 
هوا به نوعی مقدم بر شروع اثرات خاص» یعنی دیگر انواع دیگر 
خشکسالی است. شاخص‌های مختلفی برای نظارت و ارزیابی میزان 
خشکسالی‌های هواشناسی وجود دارد. به طور کلی در این مطالعه از 
شاخص‌های خشکسالی استاندارد شده بارش "5۳1 ( .۵1 6۶ ۷۲61666[ 
5 ) باران- تبخیروتعرق استاندارد شده ( 61 ۱۷16۵16-5200 
1 ,.0) 9۳1۳ و جریان استاندارد ( ,نت152 0جه دتاصدطلدز 
9 "510 برای پایش و تعیین سال‌های ترسالی» نرسال و 
خشکسالی استفاده شد. در این مطالعه از بسته نرم افزاری 1206 
تهیه شده در مرکز ارزیابیی خطرات طبیعی و برنامه‌ریزی اقدامات 
پیشگیرانه و موسسه احیا و مدیریت منابع آب دانشگاه فنی ملی آتن 
برای محاسبه شاخص‌های خشکسالی مورد استفاده قرار گرفت 
(2015 ,.۵1 21 مفاع1]). نرم‌افزار 2706 ابزاری کاربر پسند است که 
برای محاسبه چندین شاخص خشکسالی ارائه شده است. هدف اصلی 
در طراحی آن قابلیت آن در برآورد ان_واع مختلف خشکسالی 
(هواشناسی» هیدرولوژی» کشاورزی) است. 


اجرای مدل 

کل محدوده‌ی شبکه‌بندی مدل ۲ ۷۲0۳۲/۷۷۷[ مساحتی 
کیلومترمربع را در برمی‌گیرد. به‌ طور کلی انتخاب مقیاس سلول‌های 
شبکه به اطلاعات در دسترس, دقت مورد انتظار و همگرایی و 
واگرایی نتایج مدل بستگی دارد. براساس مروری بر ادبیات داخلی 9 
خارجی و نیز داده‌های موجود در منطقه. مقیاس ۵2۰۰ در ۵۰۰ متر 
برای مدل‌سازی انتخاب شد. بنابراین گسسته سازی و پیکسل بندی 
مدل 1۷]01«۳۲,0۷۷-۷۷۲ دشت مپاباد با ایجاد یک شبکه‌ی منظم 


1- 5202701560 ۳۲۵۵1۵102010۳ ۲ 

۲ )1۵ فصهاه مج 1 ممتاه) ۲۸۵۵101 02701260صه) 5 -2 
٩0162011007 12۶0۱۲ ۵ 6‏ -3 

۲ وهع01ط1 ما1۵ -4 
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تفاضلات محدود به تعداد ۳۳۳۸ سلول با اندازه‌ی سلول ۵۰۰ متر در 

۰ متر متشکل از ۶۳ ردیف شرقی-غربی و ۵۴ ستون در جهت 
شمال-جنوب با ۱۰۱۳ سلول فعال طراحی شد. با توجهدر مدل 
۱۷۲۵۱۲1/۷۷ مرزهای با جریان ورودی و خروجی به‌صورت 
بار هیدرولیکی متغیر (بار هیدرولیکی عمومی )» مرزهایی از محدوده 
آبخوان که همسو و منطبق با جهت کلی جریان آب زیرزمینی (مرز 
بدون جریان) که از نقشه خطوط هم‌تراز سطح آب زیرزمینی استخراج 
شد به‌صورت بار هیدرولیکی مشسخص " در نظر گرفته شد. پس از 
اجرای اولیه» مقادیر تغذیه حاصل از مدل ٩۷۷۸۲‏ به‌عنوان داده 
ورودی به مدل 10۳۳۲0۷۷-۱۷۷ اضافه شد ( ,۵1 ۶» هزه1۳ 
1 داده‌های تغذیه ماهانه در زیرحوضه‌های مدل ٩۷۷۸۲‏ در 
سطح سلول‌های مدل ۱0۳۳10۷۷-۱۷۷ تقسیم می‌شوند و 
سپس به‌عنوان مرز بالا به مدل ۷0۲۳۲0۷۷-۷۷۲ معرفی 
می‌شوند. مقادیر هدایت هیدرولیکی بر اساس لایه شناسی و 
سنگ‌شناسی منطقه موردمطالعه به مدل معرفی و به‌طور مشابه, 

مقادیر آبدهی ویژه محاسبه و مقادیر اولیه‌ی آن‌ها به مدل تعریف شد. 

مقادیر برداشت هر یک از چاه‌های بهره‌برداری بر اساس آماربرداری‌ها 
و مقادیر تغذیه برآورد شده توسط مدل ٩۷۷۸۲‏ در سطح زیرحوضه‌ها 
به‌صورت ماهانه به مدل معرفی شد. مدل ۱10۳۲0۷۷-۱۷۷ در 
ابتدای شبیه‌سازی به سطح آب زیرزمینی اولیه نیاز دارد؛ بنابراین» برای 
شبیه‌سازی جربان ماندگار (گام زمانی یک ماه)؛ از سطح آب 
اندازه‌گیری شده چاه‌های مشاهداتی در مهر سال ۱۳۸۸-۱۳۸۹ بعنوان 
شرط اولیه به مدل معرفی شد. بسته ,۷۷21 مدل -۱0۳۳۲0۷۷ 
1 فرآیندهای پمپاژ از آب‌های زیرزمینی را شبیه‌سازی می‌کند. در 
مدل اصلی ۲۷7۲- 1۷۲0۲۳۲/۷۷ از بسته‌های ۲5 و 5۳۲ به 
ترتیب برای شبیه‌سازی تبخیروتعرق از سطح آب زیرزمینی و تغذیه 
آب زیرزمینی استفاده و همچنین» جهت بررسی اندرکنش بین منابع 
آب سطحی و زیرزمینی از بسته 8*1۷ استفاده شد. برای شبیه‌سازی 
جریان حالت غیر ماندگان بار آبی اندازهگیری شده جاه‌های مقناهداتی 
در مهر ۱۳۸۸-۱۳۸۹ به‌عنوان بار آبی اولیه در ماه بعدی (آبان ۱۳۸۸) 
در نظر گرفته شدند. 


واسنجی و اعتبارسنجی مدل 

واسنجی مدل آب زیرزمینی یکی از مهم‌ترین و مشکل‌ترین 
مراحل استفاده از یک مدل شبیه‌سازی عددی آب زیرزمینی است. 
واسنجی مدل فرآیندی است که در آن پارامترهای نامشخص و با عدم 
قطعیت برای دستیابی به مقلدیر بهینه تغییر می‌پابند تابین ننایج 
خروجی مدل از قبیل سطح آب زیرزمینی ان دازه‌گیری و 


1- ه۲معع0)‎ ۲۱6۵0۵ 
2- ٩۵601816 ۳۱۵۵0 ۷ 


شبیه‌سازی‌شده حداقل اختلاف وجود داشته باشد. در این مطالعه در 
اقزا قل با ات رفقی ریاضین حل سابانه لاه چم شاب 
تعداد تکرار و معیار خطای قابل‌قبول برای خاتمه فرآیند تکرار اجرا 
شد. پارامترهای مدل که برای واسنجی می‌توانند تغییر یابند شامل 
خصوصیات هیدرولیکی آبخوان و تغذیه آب زیرزمینی است. در این 
مطالعه واسنجی مدل۱1010۷۷-۱۲۷۲ با تغییر مقادیر 

مجموعه‌ی پارامترها شامل مقادیر هدایت هیدرولیکی آبدهی ویژه. 
تغذیه و هدایت هیدرولیکی بستر رودخانه باهدف تطابق مقادیر تراز 

آب زیرزمینی شبیه‌سازی شده و مشاهده‌ای انجام شد. در ابتدا به‌منظور 
درک واکنش مدل به تغییر پارامترهاء واسنجی به روش سعی و خطا 
انجام شد. سپس, از الگوریتم "۳1۵57 به‌عنوان یک برآوردگر غیرخطی 
به‌منظور دستیابی به واسنحی بهینه استفاده شد ( .61 61 10۳6707 
4 به این‌صورت که بعد از اجرای اولیه مدل, در مرحله اول 
مقادیر هدایت هیدرولیکی که از نتایج آزمون‌های پمپاژ و اطلاعات 
تهیه‌شده از آب منطقه‌ای به مدل معرفی‌شده بود در حالت ماندگار؟ 
واسنجی و سپس شرایط اولیه برای اجرای بعدی مدل مها شد. در 
مرحله واسنجی و اعتبارسنجی آنالیز حساسیت هرکدام از اين پارامترها 
موردبررسی قرار گرفت. چراکه در مرحله واسنجی تغییرات پارامترهای 
حساس بسیار مهم می‌باشند. به‌طورکلی هرکدام از این پارامترها 
حساس می‌باشند؛ اما درجه حساسیت هدایت هیدرولیکی» آبدهی ویژه 
و تغذیه نسبت به سایر پارامترها بیشتر بود؛ بنابراین» تعیین و اعمال 
مقادیر هدایت هیدرولیکی» آبدهی ویژه و تغذیه عامل مهم دیگری 
برای دستیابی به یک واسنجی موفق بود. البته ازآنجاکه این مقادیر 
باید بر اساس خصوصیات منطقه تغییر کنند؛ بنابراین مهم بود مقادیر 
مطلوب در محدوده واقع‌بینانه (۲/۵ تا ۱۳ متر در روز برای هدایت 
هیدرولیکی و ۰/۴ تا ۲/۶ درصد برای آبدهی ویژه) تعیین شوند که 


4 راطعصمن۵۱ 1۱۵۷ عاماممنماعتک هه وتعع‌صلع‌ظ۲). 

در گام بعدی» پس از تخمین مقادیر اولیه ضرایب هیدرودینامیکی 
برای کاهش اختلاف تراز آب شبیه‌سازی شده و مشاهده‌ای» مدل 
تا خلت غیماتتار نا اشسفاده ادها ی یافسارم ی 
آبدهی چاه‌های بهره‌برداری, تغذیه و تراز سطح ایستابی چاه‌های 
مشاهده‌ای اجرا و مقادیر آبدهی ویژه (,45 ضریب هدایت هیدرولیکی 
(ک1)» هدایت هیدرولیکی بستر رودخانه و شاخه‌های جریان و تغذیه آب 
زیرزمینی به‌صورت خودکار و به کمک روش ۳15۲در نرم‌افزار 
5 واسنجی شد. درنهایت مدل ۷/01([71.01۷-۱۷۷۷۲ بر اساس 
مقادیر سطح آب زیرزمینی ارزیابی و واسنجی گردید. معیار خطا به نوع 


طمتقصصتای م۳2۲2 -3 
عاقاه 962017 -4 


رجا و همکاران. شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی... ۰ ۳۴۳ 


مقدار واسنجی بستگی دارد. اگر واسنجی بار هیدرولیکی باشد» میزان 
خفا مایب با مات مطالمه معطور ی شوه تیه اکتر ار 
واسنجی» جریان هیدرولیکی یا یلان آب باشده خطا جزتی قابلقبول 
خواهد بود. همچنین» مدل خطای بیلان جرمی را توسط درصد 
اختلاف" نشان می‌دهند و هر چه این مقدار به صفر نزدیک‌تر باشد. 
باشم ازاشنه تمط خقل فا قیز رات رای مطالبه بان 
ملل‌ماژی خرمادگار نخان موابد نک ووزنی ماه ارشال ان 
۱۳۸-۹ تا ۱۳۹۷-۹۸) در نظر گرفته شد. واسنجی و اعتبارسنجی 
فدل یگیب بران خومهان سور سال ۱۳۵ ۲ شچریورسال ۱۳8۹۴ 
و از مهر سال ۱۳۹۴ تا شهریور سال ۱۳۹۸ انجام و دوره‌ی تنش و 
گام زمانی ماهانه و واحد زمان روزانه انتخاب گردید. ازآنجاکه 
اندازه‌گیری تراز سطح ایستابی در چاه‌های مشاهده‌ای به‌صورت ماهانه 
انجام می‌شود گام زمانی ماهانه در نظر گرفته شد. 


معیارهای ارزیایی مدل 

برای ارزیابی کارایی مدل واسنجی شده» می‌توان میزان اختلاف 
موجود در بیلان تا برآورده شده در فایل خروجی ر مشاهده کرد 
این مقدار اختلاف باید تا حد امکان بسیار کوچک (کمتر از یک 
درصد) باشد. همچنین» از معیارهای مختلفی شامل ضریب تبیین (*) 
و ضریب نش- ساتکیف (5۳) و ریشه میانگین مربعات خطا 
(۳۹۳) و ۳8۲۸5 برای ارزیابی اثربخشی مدل و قابلیت آن در 
پیش‌بینی‌های دوره‌های واسنحی و اعتبارسنجی استفاده گردید: 


)۲ ۴( 0 - ز۳)( 0 - :0) و2 ] 2 
)۳0 كِ :۳) بجانر > )0 و0 ویر 
(۲ (:9 - )وزیا 0 
۳ ِ« ن0) رت 
[ رتش ۳ 
ر ك«ِ_ ۳ ۱ 
(۵) 100 تثِ_ِ_ِ 2۱۸5 


در اين روابط مقدار 2 تعداد مشاهدات» :0 و :0 به ترتیب مقادیر 
تراز سطح ایستابی اندازه‌گیری شده و برآورد شده (متر) و 8 و 
میانگین مقادیر اندازه گیری شده 9 برآورد شده تراز سطح ایستابی 


(متر) است. 


تحلیل حساسیت پارامترهای مدل 100۲01۷-۱۳۲۷ 
پارامترهای متعددی در شبیه‌سازی جریان آب زیرزمینی دخیل 


۳۵۲۵۵۳ - 1 
۵۵ ۷۵۵۲ 10 رمم1507009) م1" -2 


می‌باشند که لازم است میزان حساسیت هریک از این پارامترها در 
مدل بررسی شود. اگر حساسیت مدل به تغییر مقادیر پارامترهای بر 
آبخوان زیاد باشد نیاز مقادیر دقیق‌تر آن پارامتر وجود دارد 
(2011 وهمتطعل‌محطعطل؟ مهد عاوحتصتته؟1). تجزیه‌وتحلیل 
حساسیت در خصوص پارامترهای اصلی آب زیرزمینی نشان داد 
هدایت هیدرولیکی حساس‌ترین پارامتر مربوط به شبیه‌سازی تراز آب 
زیرزمینی بود (شکل ۵). میزان هدایت هیدرولیکی بر تمامی محدوده 
مدل تاثیرگذار است. همچنین» جهت و میزان جریان به این پارامتر 
وابستگی زیادی دارد. از طرف دیگر افزايش میزان هدایت هیدرولیکی 
موجب تغییر در توزیع مکانی اثر ورود و خروج آب از طربق منابع و 
مصارف مختلف در تمامی سلول‌های مدل می‌گردد. طاهری تیزرو و 
کمالی (2019 ,2۳2 0حد 11270 تعطد]) و سادات همراز و 
همکاران (2016 ,.۵1 2 ۲۲222 92021) نیز نشان دادند مدل 
نسبت به تغییرات هدایت هیدرولیکی دارای حساسیت بالا بوده و با 
آبدهی ویژه رفتار مشابهی داشته است. کوهستانی و همکاران 
(2013 ,اه 21 تصداهع00) نیز نشان دادند اگر مقادیر هدایت 
هیدرولیکی و آبدهی ویژه در مراحل تهیه مدل مفهومی به‌اندازه قابل 
قبولی از گزارش‌ها و کارشناسان سازمان آب تهیسه و معرفی شوند 
میزان خطای ایجاد شده در برآورد سطح آب زیرزمینی در اثر تغییرات 
این پارامترها بسیار پایین خواهد بود. پارامتر تغذیه نیز یکی از عوامل 
مقر در فسات سح اسغان ابیت اققا ق به‌طوری که تیه 
نسبتاً بالایی را همانند هدایت هیدرولیکی داشت. عابدی کوهپایی و 
گلابجیان (2015 مصفتط‌طهام0) 0ص ندمتامک 0ع۸0) نیز 
حساسیت مدل را نسبت به پارامتر تغذیه بسیار بالا و در حد 
پارامترهای هدایت هیدرولیکی و آبدهی ویژه گزارش کردند. 
پارامترهای که تاکنون موردبررسی قرار گرفت کل مدل را تحت تاثیر 
قرار می‌دهند و موجب تغییر سطح آب در تمامی پیزومترها در مدل 
می‌شوند؛ اما برخلاف سایر پارامترهاه ضریب هدایت کف رودخانه بر 
تعداد محدودی پیزومتر (پیزومترهای نزدیک به رودخانه) تأثیر 
می‌گذارد. این امر باعث شد میزان حساسیت مدل نسبت به تغییرات 
این پارامتر در مقایسه با سایر پارامترها کمتر باشد. 


واسنجی هدایت هیدرولیکی و آبدهی ویژه 

در مدل ۷۷۲- 1۷101۳10۷۷ پارامترهای اولیه واسنجی 
هدایت هیدرولیکی و ضریب ذخیره آبخوان بودند. اولین مرحله از 
انجام عملیات واسنجی مربوط به هدایت هیدرولیکی بود. دقت در 
تعیین اين پارمتر تأثیر زیادی در دقت سایر عوامل دارد. نتایج حاصل 
از واسنجی ضریب هدایت هیدرولیکی در شکل ۶2۸ نشان داده‌شده 


آننتت: 


۴ نشریه آب و خاک. جلد ۳۶. شماره . فروردین - اردیبهشت ۱۴۰۱ 


ات [سوبت 
هدایت بستر رودخانه تغذیه آبدهی ویژه هدایت هیدرولي 
بستر رو هی ویژ رو 


درصد تغییرات 
01۵8 0 ۳۵۲۵۵۵۵۵۵6 


میانگین مربعات خطا (متر) 
(ه) 1۱151۲ 


شکل ۵- نتایج آنالیز حساسیت و میزان خطا نسبت به تغییرات پارامترهای ورودی مدل در سه سطح ۱۰ ۲۰ و ۲۰ درصد 
۰ 270 20 ,10 0۶ ۱6۲۵۱5 6۳۵۵ ۵ )مه ون انا مدز (06مص ۵ عع‌عصل م6 ۲266 ۵0 200 فلو امه زا تاتعصهه ۵ عالیعع -5 ۲12۱۷۲۵ 


نتایج واسنجی نشان داد بیش‌ترین میزان هدایت هیدرولیکی در 
نواحی جنوبی به سمت نواحی مرکزی دشت و نیز در نواحی شمال 
شرقی (۱۱/۰-۱۲/۵) است و در حاشیه شمالی و در نزدیکی خروجی 
دشت به کمترین میزان خود (۲/۵۰-۳/۶۸) متر در روز می‌رسد. دلایل 
متعددی می‌تواند در تغییرات موجود این پارامتر در کل منطقه نقش 
داشته باشد. با توجه به این که در نواحی شمالی دشت بافت خاک غالبا 
سنگین بوده (درصد رس بالا) لذا پیشروی از قسمت جنوب به سمت 
شمال آبخوان از مقادیر هدایت هیدرولیکی کاسته می‌شود ( ۱۷/۵/6۲ 
باصعصصح۵۱0 ۱۱۵7 ماصاممتهاوتاه . مه فععصتعمص. ممتالنمومن 
4 دومین پارامتر واسنجی شده. آبدهی ویژه در شرایط غیر 
ماندگار بود. مقادیر واسنجی شده آبدهی ویژه در سطح آبخوان در 
شکل ۶0 نشان داده‌شده است. نتایج نشان داد بالاترین میزان ضریب 
ذخیره در نقاط جنوبی دشت و پایین‌ترین میزان آبدهی ویژه در نقاط 
خروجی و شمالی دشت برآورد شده است که با مقادیر بیان‌شده در 
گزارش بیلان منابع و مصارف مطابقت دارد ( عطتااتاعجمن ۱۷۷۵/6۲ 
4 اصعصمما۱۵۲۵ . ماطفصتهاونه .. فصد .. فتعع‌صام‌ه۲/ 
قسمت‌های چنوبی آبخوان به دلیل وجود تشکیلات دانه‌درشت بیشتر 
آبدهی ویژه بیشتری نسبت به دیگر مناطق دارند. چنین وضعیتی در 
قسمت ورودی جریان از سمت شمال غربی دشت دیده می‌شود. در 
نواحی مرکزی و شمالی آبخوان که لایه‌های ضخیم‌تری از رس و 
تشکیلات دانه‌ریز وجود دارد مقدار ضریب ذخیره کمتری (۰/۴-۱/۰۴ 


در این مطالعه از مقادیر تغذیه حاصل از مدل ٩۷۷۸۲‏ در 
مدل‌سازی آب زیرزمینی با استفاده از مدل ۲۷۷۲- ۱0۳۲/0۱۷۷ 
استفاده شد. برآورد میزان تغذیه به آب زیرزمینی به هر روشی با عدم 
قطعیت روبه‌رو بوده و لزوماً این مقادیر همانند سایر روش‌ها ممکن 
است از دقت بالایی برای منطقه مورد مطالعه برخوردار نباشد و 
به‌منظور تعیین دقیق‌تر آن‌ها نیازمند واسنجی مدل است. مقادیر تغذیه 
با توزیع زمانی و مکانی مناسب در سطح دشت به مدل معرفی شد 
به‌طوری که بررسی نتایج حاصل از مدل ٩۷۷۸۲‏ نشان داد مقادیر 
تغذیه به آب زیرزمینی در فصول پاییز و زمستان حداقل بوده اما در 
فصول بهار و تابستان با افزایش فعالیت‌های کشاورزی میزان تغذیه به 
آب زیرزمینی افزايش می‌یابد. به‌طوری‌که نتایج حاصل از واسنجی 
نشان داد به‌طور متوسط ۱۱/۲ ۱۸/۹ ۲۴/۷ و ۴۵/۲ تغذیه به 
آبخوان ناشی از نفوذ از بارندگی و جریانات برگشتی آب کشاورزی به 
در این مطالعه نتایج نشان داد به‌طور متوسط حدود ٩‏ درصد از آب 
بارندگی و ۲۶ درصد از مصارف آب در بخش کشاورزی به آبخوان 
| عاصمصتفاعیی مجح قععصتعص . مصلنمجمی .۱۷۷2۲ 
4 ,۵۱00۵۳060 ۱(2۵۷). 
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شکل *- ) مقادیر هدایت هیدرولیکی و ) آبدهی ویژه پس از واسنجی 
نی ماه ماه ۲۱۵۱0 »نممصرو رها) مضه رای عتاسصل وط ره) 6 6ساع۲1 


واسنجی سطح آب زیرزمینی 

انم سای حویان مانکا جرای شا فیرسال ای 6ات 
۸ موثر واقع شد. در حقیقت. تراز آب زیرزمینی در شرایط ماندگار 
به دلیل حداقل بهره‌برداری از آبخوان دارای افت قابل‌توجهی نبود. در 
چنین شرایطی» فرض ماندگاری جریان آب زیرزمینی خطای 
قابل‌ملاحظه‌ای را تولید نمی‌کند. افزایش پمپاژ نه‌تنها در مقیاس 
فصلی بلکه در یک دوره طولانی‌مدت سبب افت تراز آب زیرزمینی در 
ماه‌های فروردین تا شهریورماه می‌شود. ازاین‌رو» مدل‌سازی در شرایط 
غیر ماندگار به دلیل درنظرگرفتن پویایی سامانه آب زیرزمینی از قبیل 
تغییرات فصلی و اثرات خشک‌سالی/ترسالی نسبت به شرایط ماندگار 
مناسب‌تر است. به همین دلیل برای در نظر گرفتن تغبیرات تراز آب 
زیرزمینی» یک دوره شش و چهارساله به ترتیب برای واسنجی و 
اعتبارسنجی مدل در نظر گرفته شد. نتایج حاصل از نوسانات سطح 
آب زیرزمینی شبیه‌سازی شده با مقادیر اندازه‌گیری شده در ۲۲ 
پیزومتر در شکل ۷ نشان داده‌شده است. 

نتایج به‌دست‌آمده از هر دو دوره واسنجی و اعتبارسنجی نشان داد 
مطابقت خوبی بین تراز آب زیرزمینی شبیه‌سازی شده و مشاهده‌ای 
وجود دارد. به‌طوری که مدل روند تغییرات تراز آب زیرزمینی در دشست 
را به‌صورت رضایت‌بخشی شبیه‌سازی کرده و توانسته است تغییرات 
تراز آب زیرزمینی در ماه‌های مختلف را با توجه به تغییر میزان تغذیه 
و تخلیه از آبخوان در طی دوره‌ها و مناطق مختلف با در نظر گرفتن 
کاربری اراضی و الگوی آبیاری به‌طور مناسبی برآورد نماید. 
به طورکلی بررسی هیدروگراف‌های تمام پیزومترها نشان می‌دهد مدل 
توانسته روند فصلی و بلندمدت سطح آب زیرزمینی را در طول دوره 
شبیه‌سازی به‌درستی به تصویر بکشد. به‌طوری که اختلاف قابل‌توجه 


نوسانات سطح آب زیرزمینی در اکثر پیزومتره ا بین مقادیر 
شبیه‌سازی‌شده و اندازه‌گیری شده وجود ندارد. البنه در برخی از 
هیدروگراف‌هاء سطح آب زیرزمینی شبیه‌سازی‌شده بالا یا پایین‌تر از 
سطح آب‌های زیرزمینی اندازه‌گیری شده اتفاق می‌افتد. اختلافات و 
تناقضات موجود به دلایل مختلفی بستگی دارد. دلایل اصلی این 
اختلافات را می‌توان به خطا در ثبت تخلیه چاه‌ها همراه با چاه‌های 
غیرمجاز تغییرات ناگهانی سطح آب زیر زمینی» خطای انسانی در ثبت 
تراز آب زیر زمینی چاه‌های مشاهداتی در پیزومترها و خطای 
اجتناب‌ناپ ذیر در تخمین ضرایب هیدرودینامیک مانند هدایت 
هیدرولیکی نسبت داد. 

بررسی پیزومترها نشان داد که نوسانات سطح آب زیرزمینی 
بیشتری در پیزومترهای اراضی قره‌خان,» اراضی قزل گپی, اراضی 
گردگروی قزل‌گپی-فره‌قشلاق, گبازله و لج واقع در مرکز و شمال 
شرقی آبخوان نسبت به سار پیزومترها اتفاق می‌فتد که ممکن است 
به‌دلیل تغذیه فصلی ناشی از آبیاری بی‌رویه کشاورزان و نفوذ عمقی 
خاضل از ان تافتا طورش که پررسی‌ها پرانتانشی تفه کاریرمی ارآشب. 
نشان داد عمده محصولات تحت کشت در این مناطق شامل 
کشت‌های پرمصرف پیاز و یونجه است. همچنین, تغییرات ناگهانی 
سطح آب در برخی پیزومترهای نزدیک به رودخانه مهاباد 
پیزومترهای اراضی قزل گپی, اراضی اگریقاش, گردگروی و لج) در 
بای ههام سال لتاق افاده اسسق دلیل ایق مطلب را مب‌توام 
به وقوع سیلاب‌های فصلی و رهاسازی حجم آب زیادی در شبکه 
مهاباد طی یک بازه زمانی کوتاه‌مدت به خصوص در سال‌های اخیر بر 
اساس داده‌های رهاسازی شده از سد مهاباد در رودخانه مهاباد و نیز 
داده‌های اندازه‌گیری شده در ایستگاه هیدرومتری گردیعقوب در نقطه 


۶ شریه آب و خاک. جلد ۰.۲۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


خروجی دشت نسبت داد. تغییرات ناگهانی سطح آب زیرزمینی در 
برخی مناطق می‌تواند منجر به خطاهای منطقه‌ای شود به‌طوری که 
ما وان تور کامتل عرازم وش قاری گس راز 
محققین نیز بیان کرده‌اند سطح آب زیرزمینی در اثر تغییرات تغذیه با 
برداشت از آبخوان به‌طور پیوسته در نوسان است. به‌طوری که نوسانات 
فصلی سطح آب زیرزمینی در ارتباط با تفییرات فصلی تغذیه آب 
زیرزمینی ناشی از بارندگی, نفوذ عمقی آب آبیاری (تغذیه)» 
تبخیروتعرق و پمپاژ از چاه‌ها است (:2018 ,اه 6 76نظ 220 2۲ 
3 ,ر۵120عظ) . 

بررسی‌ها نشان داد پیزومترهای اراضی دریازء قره‌خان» اراضی 
کریت ره ماه قرل کین کوییه کووی وتان تست ی رو متن 
کمی را نشان می‌دهد. بررسی‌ها نشان داد الگوی کشت در این نواحی 
(شمال و شمال غربی آبخوان) عمدتاً کشت آبی گندم و جو و 
چنندرقند است. بنابراین» علت آن می‌نواند این باشد که تغذیه ناشی از 
بارندگی و نفوذ عمقی آب آبیاری به‌دست‌آمده توسط مدل ٩۷۷۸۲‏ در 
این نواحی توانسته افت ایجادشده در طول فصل برداشت از آب 
زیرزمینی را جبران کند. همچنین» بررسی‌ها نشان داد افت سطح آب 
زیرزمینی مداوم در پیزومترهای فخریقاه. گرگ‌تپه و سهراه حاجی 
خوش واقع در مناطق جنوب و جنوب غربی آبخوان اتفاق افتاده است. 
وضعیت الگوی کشت و تخلیه چاه‌ها در اين مناطق نشان می‌دهد به 
دلیل وجود کشت پرمصرف مانند سیب سایر باغات و یونجه و نیز 
تراکم بالای چاه‌های بهره‌برداری و درنتیجه مصرف بیشتر منجر به 
افت سطح آب زیرزمینی شده است. به‌طورکلی سطح آب زیرزمینی با 
شیب زیادی به‌یک‌باره در اکثر پیزومترها در سال‌های آبی ۱۳۹۳-۹۴ 
و ۱۳۹۵-۹۶ کاهش‌یافته است. دلیل این امر را می‌توان به کاهش 
شدید بارندگی در سال ۱۳۹۵ و افزايش برداشت از آب زیرزمینی برای 
تأمین آب آبیاری محصولات نسبت داد. 

لازم به ذکر است علاوه بر تعیین بهینه پارامترهای هدایت 
هیدرولیکی» ضریب ذخیره و مقادیر تغذیه به آبخوان رودخانه‌های 
فرعی (آبراهه‌ای متصل به رودخانه اصلی مهاباد) که عمدتاً ب‌صورت 
فصلی دارای جریان می‌باشند نیز در مدل‌سازی در نظر گرفته شد. 
اطلاعات ورودی این آبراهه‌ها از مدل ٩۱۷۸۲‏ استفاده شد. در نظر 
گرفتن رودخانه‌های فرعی در فرآیند مدل‌سازی منجر به اعمال تأثیر 
این آبراهه‌ها در نواحی مختلف دشت بر نوسانات سطح آب زیرزمینی 
می‌شود (:2019 ,۵ # تتدرتل2019۵ راه 21 و۱۷۲0 ). 
به‌ طوری که بررسی‌ها نشان داد خطای تراز آب زیرزمینی شبیه‌سازی 
شده در پیزومترهای هم‌جوار با آبراهه‌های طبیعی مانند کوسه کهریز, 
یل گوس ورین رام گر[ گیی‌و ای درک رین اس 
۴ متر) بود. نتایج نشان داد در پیزومترهای نزدیک به حاشیه 
رودخانه اصلی مهاباد[پیزومترهای اراضی قزل گپی» راضی اگریقاش, 
گردگروی» گاپیس, دارلک و لج) نیز میزان خطا کمی افزایش‌یافته 


است که دلیل آن را می‌توان به وقوع سیلاب در رودخانه‌ها و یا 
برعکس کاهش جریان آب در رودخانه در یک بازه زمانی کوتاه‌مدت 
نسبت داد. در مطالعه حمزه‌خانی (2015 ,تصهطاه۲122) نیز به 
تأثیر رهاسازی جریان‌های سیلابی در رودخانه‌های سیمینه و زرینه‌رود 
واقع در دشت میاندوآب بر نوسانات سطح آب زیرزمینی اشاره شده 


است. 

علاوه‌براین» واقع‌بینانه است که خصوصیات آبخوان در یک سلول 
شبکه ۵۰۰ متری منحصربه‌فرد است. این امر می‌تواند باعث تغییرات 
قابل توجهی در شبیه‌سازی نوسانات سطح آب زیرزمینی در بازه‌های 
زمانی خاص در مدل‌سازی شود. همچنین» بررسی‌ها نشان داده است 
اگر چاه مشاهده‌ای در لبه سلول شبکه‌بندی واقع شود تراز سطح آب 
زیرزمینی تخمین زده‌شده در مرکز هر سلول احتمالاً با مقدار 
مشاهداتی متفاوت خواهد بود (2019 ,۲167 220 21 ۷۷). 

نتایج اعتبارسنجی مدل نیز مشابه نتایج واسنجی بود که 
نشان‌دهنده سازگاری شبیه‌سازی مدل در دوره‌های واسنجی و 
اعتبارسنجی است. به‌طورکلی در نظر گرفتن مقادیر تغذیه حاصل از 
مدل 8۱۷۸۲ به‌عنوان یکی از اصلی‌ترین و مهم‌ترین پارامترهایی آب 
زیرزمینی است که تعیین آن توسط مدل طی فرآیند واسنجی. نشانگر 
میزان قابل‌اعتماد بودن روش در نظر گرفته‌شده برای تعیین میزان 
تغذیه آبخوان است. این امر اهمیت انتخاب روش مناسب برآورد تغذیه 
را نشان می‌دهد (:2016 محعت20۳0۲0 4صه 1(0۷۷/۱۸)۵0201 


2019 .]۵ #۶ ۸۱۱۷۵2۲ :2015 متتلف۲2 0جه ط212206منامکز 
20194 ,16۵0 0ص تمصع 20195 .اه 6۶ ا0م2202 ) 


خطای مدل‌سازی 

تلاش‌های زیادی در ساخت مدل آب زیرزمینی برای نزدیک 
کردن آن به شرایط واقعی انجام شد؛ اما عملا واردکردن تمام جزئیات 
موّثر بر جریان آب زیرزمینی در شرایط واقعی به مدل امکان‌پذیر نبوده 
و ممکن است نتایج حاصله تطابق کاملی با مقادیر مشاهداتی در 
محیط طبیعی برخوردار نباشند. مقادیر اندازه‌گیری و شبیه‌سازی شده 
سطح ایستابی به همراه آماره‌های عملکرد مدل برای دوره‌های 
مختلف واسنجی و اعتبارسنجی در شکل ۸ نشان شده است. 

ضریب تبیین (**) بین مقادیر اندازه‌گیری شده و محاسباتی نشان 
می‌دهد که این دو مقدار از همبستگی بسیار بالائی برخوردارند. نتایج 
معیار 1801513 با توجه نوسانات سطح آب زیرزمینی هر یک از 
چاه‌های پیزومتری نشان داد مقدار خطای ۰/۳۵ و ۰/۲۴ متر به ترتیب 
برای دوره واسنجی و اعتبارسنجی به‌دست آمد. همان‌طور که ملاحظه 
می‌گردد مقدار خطا کمتر از یک متر است که نتیحه رضایت‌بخش و 
مناسبی برای منطقه است. بررسی‌ها نشان داد ۱۶ چاه از ۳۲ چاه 
مشاهده‌ای دارای خطای ۳ عکمتر از ۰/۴ متر بودند که نشان از 
دقت بالای واسنجی مدل باشد. 


رجا و همکاران. شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی... ‏ ۴۷ 
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شکل ۷- مقایسه مقادیر شبیه‌سازی و مساهده‌شده تراز سطح ایستابی در دوره واسنجی (مهر ۱۳۸۸ تا نمهریور ۱۳۹۶) و دوره اعتبارسنجی (مهر 
۶ تا شسهریور ۱۳۹۸) برای پیزومترها در آبخوان مهاباد 
0 2009 066۵۵۵۲) ۵۵0 ماه انا صز ۱6۲۵ مااهه ماه ۵۲ مما۷۵ ۲۵0مووان 0ج 0عاهامصنه ۵ هعومجم -7 ۲121۲۵ 
منود ۱۷۲۵۳۵020 رز ممامجص۵1620 ۲۵۲ (2017 نامع ۵ 2015 ههام0) ۵100 مت)عتاه۲۷ مد (2015 هام5 


ازآنجاکه هرچقدر مقدار ۱۷۳ به یک نزدیک باشد. مدل کاراتر 
است. لذا بر اساس این آماره می‌توان گفت مدل به‌خوبی واسنجی شده 
و ورودی‌های اولیه مدل نسبتاً خوب تخمین زده‌شده‌اند. به‌طورکلی 
نتایج موّید آن است که مدل به‌خوبی واسنجی شده و قادر بوده در 
دوره زمانی خارج از دوره واسنجی یعنی دوره اعتبارسنجی نیز سطح 
آب زیرزمینی را به خوبی شبیه‌سازی نماید. به طورکلی در منطقه دشت 
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تراز آب زیرزمینی شبیه‌سازی شده (متر) 
(6۱60) عاداج1 مادنا معادانصند 
"۳ ۳ 
و و و 


۲ 
35 


اعتبارسنجی ۰ 1275 
«متاع1نا۱۷ : 


1270 
1310 1305 1300 1295 1290 1285 1280 1275 1270 
تراز آب زبرزمینی اندازه گیری شده (متر) 
(۱۵۲6۲) عاحاه؟ ۷۷/۵۲۵۲ 05۳۷6۵0 


مهاباد در اکثر مناطق سطح آب زیرزمینی کم‌عمق است. چند پیزومتر 
در مناطق غربی دشت نوسانات تراز اب زیرزمینی در حد ۲ الی ۴ متر 
در درازمدت نسبت به سایر پیزومترها بیشتر بود. این امر باعث شد 
دامنه خطا برای هر پیزومتر در فرآیند واسنجی کمتر در نظر گرفته 
شود که منجر به فرآیند واسنجی مناسب برای منطقه شد. 
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تراز آب زیرزمینی شییه‌سازی شده (متر) 
(۸۵۸6۱) 1۵16 ۷۵۱۵۲ 5۱2۱60 


واسنجی 
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تراز آب زیرزمینی اندازه گیری شده (متر) 
(6۲ع۲ع) عاصاه] ۷۷۵۱۲۵۲ ۳۷۵۵عوها0 


شکل ۸- آماره‌های عملکرد مدل برای مقایسه مقادیر شبیه‌سازی و مشاهده‌شده تراز سطح ایستابی در دوره‌های واسنجی و اعتبارسنجی 
0ص صمتاهطانای ما ماصاه) مه ۵ ماه 0ع۲موران ره 0م) اند ۵ صموزتومصی هم ی‌تاعاهاه ۵۵0۵۵ ۱۷۲۵0۵۵ -8 عساع۲1 
۵006 ۲۷۵۱102100 


۸ شریه آب و خاک. جلد ۰.۲۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


جدول ۱- بیلان آب زیرزمینی آبخوان مهاباد (میلیون متر مکعب در سال) 
۲۵۵ ۷۲۵۱۲۲ ۵2۱266 ما0۵ ت0ع و2 ۱۷۲2۳۵090 -1 م1201 


سال‌تر سال نرمال سال خشک 
۵ ۷۷۵ ۵۲ ۱۱00۵( ۵۲ 11۳۷ 
(1397-98) (1396-97) (1395-96) 


خروجی ورودی خروجی ورودی خروجی ورودی 
۱۹ 


232 326 2915 95 33-12 63 


34 2025 14 26 958 0-15 


45 0 685 0 1625 0 


مولفه بیلان 


احعدممصمع م0ص212ظ 


جریان آب زیرزمینی از مرزها 
معط ممی ما ماه صتامین 
تغذیه یا زهکشی از بستر رودخانه و آبراهه‌ها 
فتاه 24 و0ع۷۵۲۵ 1 ۵۶ معمصتدنل ۵۲ معتهطمم] 
تغذیه ناشی از نفوذ آب بارندگی و جریانات برگشتی کشاورزی 


5 ناام۲ نموه همه المگصته: ۶ ممتامماکصا ما متا معتقطه‌م۲ 


0 +35 0 2225 0 55.20 
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بیلان آب زیرزمینی 

در این مطالعه استخراج نتایج مولفه‌های بیلان آب در سه سال 
آبی خشک. نرمال و مرطوب مدنظر بود. بررسی شاخص‌های 
خشک‌سالی هواشناسی و هیدرولوژی از قبیل ۳۲ 2۳8۲ 51 
نشان داد سه سال آبی ۱۳۹۵-۹۶ ۱۳۹۶-۹۷ و ۱۳۹۷-۹۸ به ترتیب 
معرف سال‌های آبی خشک. نرمال و تر در طول دوره شبیه‌سازی 
می‌باشند. همچنین» بررسی‌ها نشان می‌دهد مقادیر بارندگی برای سه 
سال آبی مذکور به ترتیب ۱۸۲ ۳۹۰ و ۵۰۶ میلی‌متر است که صحت 
تعیین سال‌های معرف را تأیبد می‌کند. اجزای اصلی بیلان آب 
زیرزمینی آبخوان مهاباد برای سه سال آبی ۱۳۹۵-۹۶ (سال خشک)» 
۱۳۹۶-۷ (سال نرمال) و ۱۳۹۷-۹۸ (سال‌تر) در جدول ۱ ارائه‌ شده 
انش به‌ظو ‏ کلی همات‌طور گ یرای هرگنام از شال‌ها ماخ 
میزان تغذیه ناشی از نفوذ آب بارندگی و تلفات نفوذ عمقی آب آبیاری 
در هر سه سال متفاوت و قابل‌توجه است. برداشت آب از طریق 
چاه‌های بهره‌برداری» زهکشی رودخانه از آبخوان و تبخیروتصرق 
مهم‌ترین مولفه خروجی از آبخوان می‌باشند. از طرف دیگر عمده 
منابع تغذیه آب زیرزمینی ناشی از نفوذ آب بارندگی در ماه‌های سرد 
سال (دی‌ماه تا اسفند) و بارندگی‌های بهاره (فروردین و اردیبهشت) و 
نیز جریانات برگشتی آب کشاورزی و جریانات ورودی عمدتا از 


۴ ۳۲۲۵۵10۵102110 2۲011960 صه) و -1 
0۲۲۵801۲۵010۱ ۳۲۱۷۵0۵ ۳۶6۵101]2108 02701260طه) ٩‏ -2 
۴ 1۱۲۵۱۵۵ 10۷۷صوعتا -3 


برداشت آب از طریق چاه‌های بهره‌برداری 
عا۷۵ جمنلههمن طعبامعطا ممتامه‌اوماه تمه ۱۷۷ 
تبخیر از ابخوان 
نامه معط مم ممتا2۵0۵1۵ ۲۱۷ 
مجموع 
1021 


میزان تغییرات بین مقادیر ورودی و خروجی 


عمباله۷ آبامتان هه آبامرطا م۵۳6۵ مقعصعطام ۵۶ ۲2۸0 م1۳ 


مرزهای جنوب. جنوب شرقی و شرق آبخوان است. بررسی بیلان آب 
زیرزمینی نشان می‌دهد که فعل‌وانفعالات بین منابع آب سطحی و 
زیرزمینی بسیار زباد است. به‌طورکلی نتایج بیلان حاصل از مدل 
مطابقت خوبی در تعیین مولفه‌های اصلی بیلان مطابق با گزارش 
بیلان آب منطقه‌ای آذربایجان غربی برای دشت مهاباد را دارد 
| عاطمصتفاعیی هه قععصتع . مصتلنعجم . هاه۱۷۷ 
4 ,60100۳1671 1(6۷7). به‌طورکلی بر اساس نتایج بیلان» میزان 
تغییرات ذخیره آبخوان در سال خشک و نرمال در جهت کاهش ذخایر 
آبخوان» در سال‌تر در جهت افزايیش ذخایر آبخوان را نشان می‌دهد که 
این نتایج با هیدرو گراف واحد آبخوان نیز مطابقت دارد. به‌طوریکه 
هیدروگراف واحد آبخوان آب زیرزمینی نشان می‌دهد در سال‌های 
خشک و نرمال تراز سطح آب زیرزمینی به ترتییب ۰/۱۴ و ۰/۱۱ متر 
افت پیداکرده است؛ اما در سال‌تر تراز آب زیرزمینی به ترتیب ۱/۴۷ 
متر افزایش‌یافته است که با میزان تغبیرات ذخیره آبخوان نیز برای 
سال‌های خشک (۱/۶۲- میلیون مترمکعب). نرمال (۱/۴۶- میلیون 
مترمکسب) و تر (۱۰/۸۷ میلیون مترمکسب) هی‌خوانی درد 


نتیجه گیری 


برنامه‌ریزی و مدیریت منابع آب مستلزم شناخت جامع از وضعیت 
پیچیده‌ی یک سامانه هیدرولوژیکی است. ابزاری که بتواند بدون نگاه 
کلی و اجمالی به یک سامانه هیدرولوژیک بسیار پیچیده ارتباط بین 
فرآیندهای مختلف و چگونگی مدیریت و کنترل آن‌ها را در 
مقیاس‌های زمانی و مکانی مختلف مطابق با شرایط سامانه و نیازهای 


رجا و همکاران. شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی... ‏ ۴۹ 


آن بررسی کند خواهد توانست زیرساخت‌های لازم را برای حل و 
مدیریت چالش‌های موجود در مدیریت منابع آب را فراهم کند. مقادیر 
تغذیه ناشی از نفوذ بارندگی و جریانات برگشتی آب آبیاری یکی از 
مهم‌ترین مولفه‌های ورودی مدل‌های آب زیرزمینی در مناطق خشک 
و نیمه‌خشک است. بررسی‌ها نشان داده است که عدم قطعیت در 
تعیین این پارامت به میزان قابل‌توجهی در نتایج شبیه‌سازی جریان 
آب زیرزمینی تأثیر خواهد داشت. استفاده از مقادیر تغذیه با در نظر 
گرفتن تفاوت برداشت از آبخوان و آب سطحی برای تأمین آب آبیاری 
در سطح کاربری اراضی در مناطق مختلف و درنتیجه نفوذ عمقی 
به‌جای لحاظ تغذیه به‌صورت درصد یکنواختی از بارزندگی و کل آب 
آبیاری در سطح آبخوان برای تبیین و تجزیه‌وتحلیل شبیه‌سازی روند 
نوسانات سطح آب زیرزمینی در سطح آبخوان موردتوجه قرار گرفت؛ 
لذا در این مطالعه در مدل‌سازی آب زیرزمینی آبخوان مهاباد با استفاده 
از مسدل ۷0]«۳۲,0۷۷-۱۷۷۲ از مقادیر تغذیه حاصل از 
مدل ٩۷۷۵۸۲‏ استفاده شد. 

نتایج تجزیه‌وتحلیل حساسیت در خصوص پارامترهای اصلی آب 
زیرزمینی نشان داد هدایت هیدرولیکی حساس‌ترین پارامتر مربوط به 
پیش‌بینی تراز آب زیرزمینی بود. پارامتر تغذیه نیز یکی از عوامل موّثر 
در نوسانات سطح ایستابی است. به‌طوری‌که حساسیت نسبتا بالایی را 
همانند ضریب هدایت هیدرولیکی داشت. به‌طورکلی برآورد مقادیر 
تغذیه آب زیرزمینی به هر روشی با عدم قطعیت روبه‌رو بوده و 
به‌منظور تعیین دقیق‌تر آن‌ها نیازمند واسنجی مدل است. بررسی نتایج 
حاصل از مدل ٩۱۷۸۲‏ نشان داد مقادیر تغذیه به آب زیرزمینی در 
فصول پاییز و زمستان حداقل بوده» اما در فصول بهار و تابستان با 
افزايش فعالیت‌های کشاورزی میزان تغذیه به آب زیرزمینی افزایش 
می‌یابد. نتایج حاصل از واسنجی نشان داد به‌طور متوسط 7۱۱/۲ 
۹ ۲۴/۷ و ۴۵/۲ تغذیه به آبخوان ناشی از نفوذ از بارندگی و 
بهار 9 تابستان اتفاق افتاده است. 

به‌طورکلی بررسی هیدروگراف‌های تمام پیزومترها نشان داد مدل 
توانسته روند فصلی و بلندمدت سطح آب زیرزمینی را با توجه به تغییر 
در میزان تغذیه و تخلیه از آبخوان در طی دوره‌های مختلف تنش در 
شرایط غیر ماندگار را به‌طور مناسبی به تصویر بکشد. به‌طوریکه 
اختلاف و خطای قابل‌توجه بین مقادیر سطح آب زیرزمینی شبیه‌سازی 
و اندازه‌گیری شده در اکثر پیزومترها وجود ندارد. بررسی‌ها نشان داد 
تغییرات ناگهانی سطح آب زیرزمینی در پیزومترهای واقع در نزدیک 
مرزهای ورودی و خروجی آب زیرزمینی یا نزدیک به رودخانه مهاباد 
می‌تواند منجر به خطاهای منطقه‌ای شود به‌طوریکه مدل نتواند 


منابع 


به‌طور کامل شرایط واقعی را شبیه‌سازی کند. علاوه براین دلایل 
تلف ارحیاف زر یم یه مها راو جا ای خیش ار 
خطای انسانی در ثبت تراز آب زیرزمینی چاه‌های مشاهداتی در 
پیزومترها و خطای اجتناب‌ناپذیر در تخمین ضرایب هیدرودینامیک 
مانند هدایت هیدرولیکی می‌تواند در بروز خطا بین مقادیر اندازه‌گیری 
و شبیه‌سازی شده نقش داشته باشد. بررسی ننایج بیلان آب زیرزمینی 
آبخوان مهاباد نیز نشان داد میزان تغییرات ذخیره آبخوان در سال 
خشک و نرمال در جهت کاهش ذخایر آبخوان. در سال‌تر در جهت 
افزایش ذخایر آبخوان بود به‌طوری‌که این نتایج با هیدروگراف واحد 
آبخوان مطابقت داشت. 

نتایج این مطالعه نشان داد استفاده از مقادیر تغذیه برآورد شده 
مدل 81۷۸ به‌عنوان یکی از اصلی‌ترین و مهم‌ترین پارامترهایی آب 
زیرزمینی می‌تواند امکان شبیه‌سازی و برآورد دقیق‌تر نوسانات سطح 
ایستابی و ببلان آب زیرزمینی توسط مدل ۱010۳1,0۷۷-۲۷۷۲ را 
با ابعاد مختلف مکانی فراهم کند. به‌طوری که نتایج مدل تلفیقی 
نشان می‌دهد مقادیر تغذیه برآورد شده با در نظر گرفتن کاربری 
اراضی و تأمین الگوی آبیاری از منابع آب سطحی و زیرزمینی در 
مناطق مختلف به‌خوبی توانسته است نوسانات سطح آب زیرزمینی را 
در موقعیت‌های مختلف آبخوان شبیه‌سازی کند. به‌طور مشخص 
بررسی‌ها نشان داد افت مداوم سطح آب‌زیرزمینی ایجاد شده در 
مناطق جنوب و جنوب غربی آبخوان (پیزومترهای فخریقاه. گرگ‌تپه 
و سه‌راه حاجی‌خوش) به‌دلیل وجود کشت پرمصرف مانند سیب و 
یونجه و نیز تراکم بالای چاه‌های بهره‌برداری و درنتیجه مصرف 
بیشستر بوده است. به طور کلی از مدل تلفیقی -۲ 5۷۷۸ 
21 101۳1/01۷ می‌توان به‌عنوان یک ابزار مدیریتی در 
تهره‌برداری از -سایع آب تلقیقی استفاده کترد اعسال ستاریوهاق 
مدیریتی مناسب و تبیین نتایج آن می‌تواند درک صحیح‌تر و 
واقع‌بینانه‌تری از شرایط کنونی برای مدیران فراهم کرده و ابزار و 
رویکردهای مطلوبی را در اختیار آن‌ها قرار دهد تا به پیاده‌سازی 
راهکارهایی بپردازند که حفظ منابع آب و درعین‌حال توسعه‌ی پایدار 


در دشت را در پی داشته باشد. 


سپاسکزاری 

مولفین از مدیریت و کارشناسان ستاد احیای دریاچه ارومیه که 
امکان دسترسی به اطلاعات موردنیاز تنظیم این مقاله را فراهم آوردند 
مراتب قدردانی و سپاس خود را اعلام می‌دارند. 
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معا ماهصمتع مه ماه مساو ۵ ممتاهاندته عتامصزطمن :2016 هو صهتل0عمص20 م4صه ریق ۲(0۳۷۱۵)۵0201 . 


:۴ .485-496 :(20)01 وطلهعمتوص ات ۵۶ امتمژ انم 0مصفماهه ۳۱۸۵02۵020 هد ۱۷۲۵۲۲۲/۵۷ 2208 5۱۷۷۸۲ 
.10۰1007/912205-015-0354-۰ 


ص ۵08 5۷/۸۵۲ حصم نع ومیتله۷ موتقطهع۲ ۵۶ ممتام۱/0112 .2015 هبو صقیلمعمعهمم هه ریق ۲۱(0۱۳۷۱۵)۵0201 . 


1 24 6اه ۷۷ .صتقام ۳۱۶۵۵2۵۵20 ۵ سمل بمهه نامع ماقانصته ‏ م۵ ۰ ۱۷۲۵۱۲۳۲۲۸۵۱۷ ۶و . اع0مص . لهع)مصصعط)ه۵ظ 
۰ 02 8869/6200 10۰1 :۱۱۵۲ .(200اوهاه طفتاعصظ 1 صهتونع۴ 12) .337-348 :(1901 5016860 

۱ 110۳۶ مممیته هم ععفننامه ۵۶ فعمصه۷ نامع عطا مامعتاوه۷ص1 ما عصتام0مهه مماه0ا0۲۵ ,2015 ۴۰ تصحطاطممصه۴۱ 
۰ .(صتهاج. طاهمصهن۷ :او عقعع) «ومامصطع مصتعصعه مامصعر مصتوتا فمتمتامم صمتام)امامه ماگنه مصترامم2 
(صهتوع۲ صل) ۵۰ 175 ررع10مصطمع1 ۶ اتویعب تلا کتتقطاه روتفعط ]1 

۶ ۸۵۱1621008 .2012 هر هتصصهللمطوه]۱۷ مه ,۱۷ نطوع0ه ر.ظ صقصه‌طهدات ره طع220له بنک وته7 1۳2 ببظ 12207 
:(2063 فتاه رعمامعع۳۱۷۳0 .موی رصتقام 0۲محاقطوتها مط صا دامع تمه نامع )عت0هعج مه عصتام0مصظ. هام-امصو۵" 
۸05://01.0۲2/10۰1007/510040-011-0814-۰ظ . 435-447 

نامه عطا ۵۶ وتورلمصه بوتاتطمصنمافداه 0۵02۵10201۷6 .2018 ٩.‏ 129284 4ص را طاتطحطتصوظ ,بط صهل0همطوما صهتمهکز 
۳۵۹122 18) .1587-1601 :(31)6 5011 200 تملج ۷۷ ۵۶ نامز .متام 0صهززه-صهوهعمطکا ایام :تاو عقهع) مصمولاه 
۰ 0 :1۲۱۵۲ (اععتادطاه طوتاعورظ 


17 


18۰ 


19۰ 


]5۷۷۵۸ 22160عاص1 عطا ۶ه ممنای‌ناممه مه اصمجصم۱۵۷۵۱0 .2008 .1.0 ۸010 20 ررو.۲ ۷۷۵۵ ,۲.۱۷ عمط ,۱۱۰۷۷۰ نک .20 


0-۰ [/۳۸۵5://01.:0۲8/10۰1016 1-16۰ :(356)1-2 بوما0ع ۴۲ ۵۶ تمصیامول .ام0مهر ۱۷۲۵1۲۲۱۲0۱۷۷ 
عصلود ۱6۷۵ هدمع ۵۶ میمصت لهتنممصتا ار :2013 هر تصهطعطهنا مه ری تطعقاقط طه۶ع۱۷ ر,.لظ تصهاعع‌طاهمک 
00/۱۵ لمصمتامصعاص1 رمع مصمافماهت رلعطمعماه ممصدا :0و مقهه ه) ۵۴و ۷۲۵۱۳۲۲/۵۷ 

858-۰ :(108 موم [0162ع10ظ قمع اهع1ع1010ظ 
ما ۲6۶۱165 حجمتامناحه تمه 268و لوا ۵۶ «مصمتله معط گم صمتاهناله۳۲۷ .2015 .۲ زبتوداط مه .۱۲.۳۲ ط2162206ظ۱۵مکز 


21۰ 


22. 


رجا و همکاران. شبیه‌سازی بیلان آب زیرزمینی با استفاده از مدل تلفیقی و جامع آب سطحی و زیرزمینی... 


۵ 238] 10165 صا 1062160 تمکنامه ب«مللقطه عمط‌لماته صا ممتامانته ‏ ماج صتامع ۵۶ پرمحتتاهعج ‏ فط ۰ 1۳2۵۲0۷۵ 
۰(امحتاوهه .. فنص ده صقتقعظ ... صّ)... ,129-138 :(1403 متام . عصتممصتعصظ .لت وعته0ه]۷ .رصهع1) 
6-865-1۰ 1 1۳/27016016 26 .12۲65 مق . ز۱)0://۴۲۵۵ظ 

۷ ۵۱6 ۷-صناهع ‏ معصعععرل-عاتصظ تمطمتعصمصصیل-ععطا. عقانالمصصر هر .1988 ۷۷۰ طمییه‌هایم۲۱ مه ,م۱۷۲۵ 210مم0ن۷]6۱ 
۰ 0۵0102102 .۱۲۲.9 تفلصع ۷۱۲ رطماوفک .6 امظ رقهمتاهعتا۵ع1۳۷6 انوم عملج/۷۷ 0۶ ومبامتصطاهع 1 ,2۵061 
۰ و( 

۶ ع۲۲۵6۵601 :1 .وملهعو مصص1 مامت اد عصتهماتصمجظ اهیناریا :1995 .1 ]19ع(16 مصح ر.ژ.اظ صملوع1(۵ ر.ظ.1 عمکزم]۷ 
.233-۰ 900160۷ اهمنع۷]۵)۵0۲010 صهمزمصصه :12 رفهالهنا .102مامصصنات 0عتاووض ج۵ عمممبع‌تومن ط90 عط 

۵ وعماعع1 منصمع000عطاصه مط) ۵۶ )دمص فطل ۶م صمتامتتلد 7 1۳ :2018 ۲ طع12029220 صه ریک تعمصصله ره 4عطه] زم]۱۷ 
۶ اهب نصا مصتهء‌صتودط همم مه 0۷۱ ۵۶ لمصتاول .019 مه عصتعصمو مامجصعع عصتودا ولفتت هتصصینا معام [ 
موه طفتاعصظ ۳ جهزوه۳ 18) .59-70 :(4802 1207127 

۵ طا ص1 عوتمط0ع۲ ماه صامم آمتصمامم عصتام۱۷۲۵0 .2019 ,1 «عاتظ مه .و یط ربا عفصعاطهتطظ ر.ظ مفووم۷ 
۰.۰ ۰ هاگ ...۰ ۵ امن ۹۷۷۸۲۰۱۷۲۵۱۲۱۲۲۸ . طا91 .. صلقفجه... هلان۴ 
.2٩0.20 19۰‏ [/16 10.10 /۱)005۹://001.0۲2 

مصتاهانله صا 0زمز ومع گ۵ عمصهاع0م1 :2011 .۲ زموننج۳ 20 ربسا.ظ مممصطمو رال ما۷۷ .۷۷ عصنک یه عنهل 
-1285 :(47)6 (1۵۷۷۴۵) مهله50عظ ۵90۷9۲۵۵ م۱۷۷2 مومتنمصصد لول فلوم مم)هلصته تفتان مه م۷۷۵2 
61 -۸10۹://01.0۲8/10.1111/1.1752ظ .1297 

7 ۱۷۲۵8۵8۵6۴0606 90۱۲۲۵۵ عماج ۲۷۷ .۲۵۷191)60 )نامب لهعتوما070 بط ۵۶ اممطرووعووظ .2009 06۰ 152715 هه ربا فتاصه‌طالهل 
1269-008-9305-۰ 05://801:078/10.10017/91 881-897 

ص کععصقطه ععه/ مصنامع )منعت صرح 2۵لمصه 1۵ وعتیعو عصصتا ۵۶ ومطامصظ عنام ۵۶ مموتنهم‌صطمن .2019 بر تبعووهآ 
۳۵۲۵۹1۵8 1(2) .58-68 :(13)01 ممهصته:نا مصج جمتامعت1 ۵ تفصتهژ صوتصوا ,رملق یتنا ۵۶ طامص) صتقام 0صفته]۱۷ 
9اعحتاوطاه طوناعط۲ 

۵۵۵۵ 602۷1۱0۲۵1 ما ومتصصمصول صموود ۵ عوتا م1 .2014 هر عوطاه)1 9212 مه ریبک مه طامصته مه رای دوه 
0 .937 :(13)04 لوه۷ ۱ مک ممممزهگ ۵۶ تمصیمل مدای بصتقام متعطاه1 گه ممتاجهمتصماع ۷1۵۱۵ ۹216 200 
اصص 1644-12 -11/2۲016016-1 .1۲:26 :6۵۰/861 انا رصهتوبه۳ 

۶ 0۵ رجمتاعناتافجمی حصصهل ۵۶ ]عدمرص فطع 160مصصفده منز ل2صمطه]۱۷ ۵۶ فتورولممه لمعومامطام:۱۷۲۵ .2015 ,۲ ۱۱2۷۷۵۲۲ 
:۱ موه طفناعص اد صهتویعظ ‏ صل) . ,155-178 :(15)37 وعمصملمو. . لهمنطاممتعهعن 0‏ ظا. وعطم‌تهعوم۴ .۰ 1160ضج خر 
6 ۵ )7 

۱۷۱۵۱۱۳۲۱۵۷۲۰ 10۲ ممتاحنم؟ وماسمه خر ,۲ ۰۵۱۷۷ ۱۷۲۵۲۱۱۳۲۸۵۱۷۷ .2011 ,۱۷۲ ۲0۵/۵1 مه ریق همه تیک بممممبس و1 
۰ ۵/0102 6-۵37۰ ۱۷۲60009 مه فعی‌متصطامع] وتا و۲79 :2005 

,206 2عع0صمصنطنسجمم! مطا 0۲ 061مصظ ما تمه نامع ملدعو-اهدمنوعع خ .2007 اما «علتصوظ 4ص .۲ ماو 
۳)]۴05://001۰0۲2/10۰1007/510040-007-0197-۰ .1457-1472 :15 متا[ زع۵0۱0ع۲۵ ۳۱۷ ,هناع۱۱16272 

2610064 زگ ۱۷۲۵06 عمتندمع:ظ ۵۲ ممتاه)صمصصینن بت ه1]۵۲1 ۸۲-۱۷۲۵۲۱۲۲/۱۷۷ ٩۱۷۷‏ .2017 .۲ «عاتوظ 0ص2 ر.و ابو 
۰ :۱179۸ ,00 رفصتلامن ۲0۲۲ :الوهب نا ماگ 240هام رعطتهعصتومط آماممصمم صظ مه ذت) ۵۶ )همحانهمم(1 
.۰ -۱)005://001.0۲2/10.1111/1752 

,ام مطتامطتاوه ما اع0مصه ]۵۷۷۲۸۵۸ ۵۶ صمتاه‌طانله موممم۷]۵1۵ :2021 ۷۲ تطفتدژه 1 هه ,۷۲ ۵0[مصتونهط ربنم وزه۴ 
8 .11-34 :1704 ط۲دعوم1 وعمسا0وم۴ عماج ۲۳۵۵-۷۷ .رصتهام 24طمطه]۱۷ :باه عفهه ع) ۷۱۵۱0 وم مضه رصمتاه ام فصهاممهم 
۰ 0 :]۱۱۵ (0مماوطاه طفناعص طاز صهتونه۲ 

صنقام 2متنطگ که تمذهطه‌ه متلنیجیل بط مق وتوروامصه با تاتقصهو ,2011 0۳ ومطمل‌ممطوطف همه ریکایه تاموتا تفگ 
(امجتدهاه افتاومط طزه معتوع۴۳ 18) .102 :(22)78 ۷۷۵256۵۷۷۵۲۵۲ 200 ۷۷2۵۵۲ 0۴ لهصتا0[ ,۲۱۷۲۷۷۲۵۲ عصتوا تعتنا20 
۷۵۵ 13 ؟۶ه عمتامعع(1 ملامع)0ظ .2016 ۱۷۲۰ تطیامصت)ومطهاط مه ر. متفگ ۷.۴ حصهل۵مطوع۱۷]۵ ۱۵2281 
0 عصللهم۱ فصتفوظ :رتاو عققت) مارگ ممتامصمصا متطامدنومعن0) صا دقعم۶۵ظ ۱6۵۲ متا امه مصلونا وععباهوم۴ 
۱ .2090200 طفتاعحظ طازه صقتقوهع۳ م) .379-389 :(303 ۳۱۷۲۵0۱02 معط ۵ لقصنمژ. . صوتصهیا برصتقاظ . مترطاه ]1 
6۰ 1022059/1۳ 

عصناع۱۷۲۵0 ۲۱۱0۷ ماه ص0۵ فصح ممه]تاو ]18 .2019 یه و122 20 ,.لانظ 2261 بخ ۸2206 ربق کنا52202100 
2۴ (1) .521-536 :(33)04 011 4ص عاه ۷۷ ۵ ماو 6۷۲۵۲۰۱۷۲۵۱۳۲۸۵۱ ۷10 0عطوتمج ۷۷ 0۲مطاقطاوزعاظ ط1 
۰ 0۷۷ :۱۱۵۲ ()26او20 طوتام‌د۴ 

2 متفه ملد مه ما ماه صنامین عصتام۱۷]۵0 2017۰ .۸۵ ۱۱201 مه ره صهل)فطوه]۱۷ تتقطعفه ریق عطمته‌طامو 
۳۵۲۹120(۰ م) .41-50 :308 صع ۵ ماو ممماهتان راملما موم و0 عملونا تعننمه صتقاظ تقاوه‌مفط 
۰ .5۷7۲۳2021 [05://:,01.0۲2/10.22059/1) 

۶ اصنهاتهعصنا منهج ۵ اصعططروو6و۵ظ .2016 ۰ ,۷۲ تهعمااظ 2۵عت۵ظ مصج ریخ تناممتجطاخ ر.ظ مها )9202 
61-۰ :(22)6 20100 00896۲۷) 6011 م4صنج عماج ۷۷ .رصتقام 4صدرتاظ :مناد عقفون) 00طامجط نانام متفه عم ,]۷۲۵ 
۰ :۲۵۲ .2090200 فاعم اه صهتوهظ م1) 

1 .۰2006 .۲ وتمصصصت. مه ۷۷۲۰ وصتظ ررظی ورهت) ونیا هلا ملظ )صنا وک فک میکا.ظ صملصهعو 
3335-0۰ :(20)15 متام لقعمتاح‌صعاما مخ .عصمتوع: تج مه 0تمتصعو طا مفتمطهع مه نامع ۵۶ ولوم‌طاطره 
۰ -0)0۹://001.0۳2/10 

ص فصمتامهتعاصا عماه ماو گ اتازصامتنه۷ له۵۲مجصه-متاهمو ۵۶ فلورولقمه .2019 ۱۷۲۰ طعمک مه .]۷۲ تنگم 
:(78)6 و۵مطمز0 اتف ماممصمم صظ باع0مصه ,۷۷۸۲-۱۷0۵1۲۲ ۱60منتاهه 2 عون ریصقت رصتفده ۲۱۷6۲ ومفطه‌عقطن) فط 
2665-019-8206-۰ 10۰1007/51 /605://001۰0۲2 201۰ 
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۲ نشریه آب و خاک. جلد ۰.۲۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


5۲ 4۶هام ممتامصتادء هدمع مماه 0و ۵۴ صمولنهم‌مصم) .2009 .5 عاممطصملطنا فصه ری).۲ تعصصماظ ر.1 42صعطان 
1427-1۰ :17 0۵1[ ۷ ۷( 7 ,21762 مه .۰.۰ همه عص 
.,-۱۲05://00۵1۰0۲2/10۰10017/510040-009-0445 

م۵ 2۲0۱۲۱0 20 0مطومزه 0مصاطصمی ۵۶ ممتاهعتامج مصه ب«وما00مط۱۷۵ .2000 .ظ. فمناعظ همه ۷ فامعهمطمم 
16/50022-1694)00(00293-۰ 10.10 /۳۸۲۵5://01۰0۴8 185-201۰ :(236)3-4 ۳۱۷۵۲۵۱۵2۷ ۵۶ اهاز .عحعصهک ما و۵0 
۰ ,1۵529 طا تعکننامه فصتقاط طعنظ مطا ما اتلزمصتهافته مه موتمهاهعر ماه تمعن :2005 ۱۷۸ فیامعاهمطم۵ 
.-۵5://001:018/10.10017/510040-004-0385ظ :351-365 :(13)2 لمصیامل «ع010عع1۲0 ۲۱۲ 

فصتحاظ 1۷۵۲۵۵01 مط) ۶ ویعه/200عظ معط صا موعقطامع: عماه7صنامتع «مااقطه عصت‌مصصتوظ ,2005 ۳.۹۰ طفتصرمت هه ری من 
۰ 10 /0۲05://001.0۲2 .795-807 :(19)3 ۲۶۵689868 ۲۱۷۲۵۱021621 5۷۷۵۲۰ ولو 

مه نامه ص1۱0 مد اع0مصص ,۷۲۵۲۲ ها عصتاممصظ بماهصتا0ع0 :2019 ,۱۷۲ هک 20 ریق ۰11220 120611 
0 .89-104 :(12)40 عطاتعمصلوطظ! عماه ۷۷ .عممتاتصمی معبفبظ لمح اصعوعيم لصا متقاو لهمنومامم1۶0۵ظ که ممتامناهبم 
139812۰ 20086377. ۰.1 20.101 :1۱۵۲ (209072600 طوناعص طز صهتوعظ 

۰ 2۱6 1۳0۱ منم 6و2 فلورلممه نمی عم مهو ه هرن :2015 .0 واتللمو۱ مص2 ,بط ولام‌عصه ۷ .نا فقعلمز 1 
.145-014-0178-۷ 10۰1007/12 /۳۲۲۵5://001,0۲82 .697-709 :(803 وعلاقصصملصا عمصمته و ط) ت۴2 

وصتعام 20طقطه۱۷ ۵۶ موب ۷۵امصنزومی تقصصتامه ۵۶ [06معظ 0۵2۵01)201۷6 :2018 ٩.‏ طقصهم‌صهل ۷2۸2 هه ر.ظ صهعم21120206 ۱۷ 
م۴۵۳1 182) .631-640 :(50)4 عصععمصلمط۴ 2۷1 ۵۶ امدربامژ ایادخ رفعمتاموع؟ میج 0صامعرع ما نج عمقتاه 
۰ 10.22060/01۳۳7 :۱۱۵۲ (260اوط2 طوتامط۴ 

1016020-010۸ و.ظ 6[602 ۱۷۱۵0142-1 ریک مامداط ریا مصمصصتت رخ 1۳۵000 ریز.ل مصمزم1۲۵۵۵2-۷ رب.5 مصحتهقعاهمع] ۲۷ 
۷ ۱6 ط1 ۲۵80۲۲6۵5 ماه م24 فاطمتامت و0 فعووعع۵ع۵ عصتمطنه/۳ ۵ ۳1۲6۵۱5 .2011 .[ 2۵02122 220 ر.گ ۱16۲12عظ ر.ژ 
۰ 0 .203-212 :48 62700و۱۵ ماقصتات .(1930-2006) هابافمتمعظ صهتنه0 

۶ ۵12266 ۲۵90۱9۲۵8 ۱۷2/۵۲ ۵۶ وعتلتتاو م۲۵02 .2014 .اصمصم۱۵۷۵۱۵ ماطمصتفاودی مه فتممصتعصظ عمجم ماه ۱۷۷ 
.0 81 دنه 00۲و 0دصمطه]۱۷ ۵۶ معصهاهها وعمعاموع1 معا۷/۵ 0۶۴ ۲0۵0۲۲ رههته اصعحصطماق معلضا متصرنا ۵۶ ففعته رنه 
(۳۵۲5121 

فصصماوره گنه داد 1۳۳۱22160 باع۷تعصماصاً صز وعومممروع؟ مدرد 0۴ )حعصصووع۸۵5 .2019 .151 ب«ماتحظ 20 ریک ۱۷۷۵۲ 
۰ 0 1033907۷7۱ /۵9://001.012 1576۰ :(11)8 2۸۵۲ ۷۷ ۰ ۰۵۱۷۷ ۵۲-۱۲0۵1۱۳۲ 5۷۷ ولو 

۲ ماه ۵۴ مملاحانصله ۷۵ تقصمطام‌م‌صصمن .2018 ۷۲۰ زماهعه 20 وت).1 ۷۷1016 ,.۴5:۷۲ ۵60۲08 ,. 1 21167ظ ریک ۷۷۵1 
۰ ۱۲۷-1 ,9۷۷۸۲-۱۵۲۳۲ 0ععاصنا . فمطا . عصلوت.. مصماووه . متعمامعظ ...مه یادهاناوعمه یه متام ظ1 
۰ 1 ۹01.20 10۰1016/0۷ /۱)0۵5://01.0۲2 1-10۰ :122 م0102 0صه عصتاام۷۵0 1ه)ممصصطمتدطرظ 

مد کصمتاعمعاص ماه م0۵ -عمه تاو لمح موتهطهع۲ ممهصناوتع رقعووعع0۵عم لهعنع۲۵۱۵ ۳۲ .2010 .۱۱.9 م۱۷۷۲ 
من رمنهع۷۷۳ .5 ۲۱۵۷۷۵۲۵ .60 186 رقههع۸ نع مه 0تنه ما عصتام0م]۷ تماج0صا۵زن .قههیه 96-2110 20 2110 
.5-۰ ,۳65 ات۱0۱۱ 9 // ۱۱۱ 11 ۹ 230 ۱/۳۹۹ .۸ 
۰ و ( 

۵۵ ۱۳۱۳ ۱۲۵۱۲90 همه ممطبمو ما ممتافتایی معتهطهع مه 0100 :2003 ۳۲۱۰۴۰ صقحصنامعظ 2020 و۷۲ 2۱۲ 
.105://001.0۲8/10۰1007/510040-016-1898-۰ .۰ 92-9220-000-3 ر۲18ظ 
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